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ABSTRACT 


The external nucleate pool-boiling heat-transfer coerfi- 
cient of a horizontal smooth copper tube in nR-114-oil 
mixtures (0 to 10 percent oil) was measured for heat fluxes 
from 1 to 100 k¥/m2 at two different saturation temperatures 
(52.2 06 .0n40 060 A A copper-nickel tube coated with the 
Union Carbide "High Flux" coating was similarly tested. The 
High Flux coating was found to improve the heat-transfer 
coefficient by at least a factor of 7 in oil-free R-114. Oil 
resulted in about a 20 percent reduction of the  heat- 
transfer coefficient of the High Flux surface at heat fluxes 
less than 30 kW/m? and up to an 80 percent reduction at heat 
fluxes above 30 kW/m? with greater than 6 percent oil. 
Under all conditions, the High Flux coated tube outperformed 


the smooth copper tube. 
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I. INTRODUCTION 
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A. BACKGROUND 


The 7 Navy currently uses refrigerant ) in 
centrifugal chilled-water air-conditioning plants atoard 
submarines and surface ships. The Navy hopes to reduce the 
Size of these units and increase their performance by using 
enhanced evaporator and condenser surfaces. in experinent 
by Arai et al. [Ref. 1] produced a prototype 200-Ton RUE 
centrifugal water. chiller that vas 28 percent shorter in 
length and had a 50-70 percent improvement in the overall 
heat-transfer coefficient by employing the enhanced surface 
"Thermoexcel E" made by the Hitachi Company. Comparisons o£ 
various enhanced commercial tubes by Yilmaz and Westwater 
İRef. 2), Marto and lepere [Ref. 3], and Carnavos [Reta 
for various refrigerants other than R-114 indicated that a 
porous-coating-enhanced surface, such as Union Carbide"s 
"High Flux," will exhibit the best Boiling heat trans ии 
performance in a pure refrigerant. 

The High Flux surface (see Figure 1.1) consists ofa 
Sintered metallic matrix bonded to a metallic substrate. 
The surface is produced by coating a smooth tube with a 
binder-solvent mixture and then applying a mixture of metal 
and braze alloy powder; the tube is placed in a furnace to 
evaporate the solvent, binder, and melt braze alloy thus 
forming a porous structure having multiple reentrant cavi- 
ties to enhance nucleation. (Ref. 5] 

Since the High Flux surface will be employed in a 
refrigeration unit using ап/"611-ТиБт иси cos hernetically- 
sealed, compressor, some amount of oil is always present in 


the evaporator. Studies by Henrici and Hesse [Ref. 6) for 
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smooth tubes and Stephan (Ref. 7] for the Gewa-T surface 
(manufactured by the Wieland Company) in &-113-0oil mixtures 
indicate that the heat-transfer coefficient of enhanced 
surfaces can be significantly altered by oil. 

Experimental data showing the effect of oil concentra- 
tion on the heat-transfer coefficient of the High Flux 
surface in R-114 are lacking, thus motivating the present 
investigation. This investigation was funded by the David 
W. Taylor Naval Ship Research and Development Center. 
Details of the experimental apparatus are described by 
Karasabun [Ref. 8]. The smooth copper tubes were supplied 
by the Wieland Company. The High Flux coated copper-nickel 


tubes were supplied by the Union Carbide Corporation. 


B. THESIS OBJECTIVES 


The objectives of this thesis are: 

1. Take boiling data on a smooth tube in R-114 with and 
without oil for comparison with the data of otier 
researcners, and to provide baseline data for evalu- 
ating the boiling performance of the High Flux tube. 

2. Take boiling data on a High Flux tube for various oil 
concentrations (0 to 10 percent by mass). 

3. Study the effect of saturation temperature on the 
R-114 boiling behavior. 

4. Attempt to sample oil locally in the near vicinity of 
a tube to investigate the possibility of an oil 
concentra tion gradient around the tube during 


operation- 
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11. REVIEW OF REFRIGERANT-OIL MIXTURE BEHAVIOR 


A. WUOCLEATE BOILING OF REFRIGERANT-OIL MIXTURES FROM SMOOTH 
TUBES 


In 1963, Stephan [Bef. 9] published a milestone paper on 
the influence of oil on the boiling heat transfer of R- 12. 
The effects he noted have been observed in most refrige- 
rants, including R-114. In 1972, Henrici and Hesse [Ref. 6] 
updated Stephan's work for R-11ü-oil mixtures boiling from a 
smooth copper tube. Figures 2.1 and 2.2 summarize Henrici 
and Hesse's results. Figure 2.1 shows that oil generally 
lowers the heat-transfer coefficient, and that at high heat 
fluxes and high oil concentrations (10 percent), the effect 
grows more pronounced (slope decreases). Figure 2.2 shows 
that at some oil/heat flux combinations, the heat-transfer 
coefficient may actually be improved by the addition of oil. 
Chongrungreong and Sauer [Ref. 10] suggest that the heat- 
transfer behavior of refrigerant-oil mixtures can be attrib- 
uted to 5 major factors: 1) the physical properties of the 
refrigerant-oil mixture, 2) the saturation temperature (or 
boiling pressure), 3) the tube diameter, 4) the surface 
condition of the tube (roughness), and 5) the hydrostatic 


liquid head above the tube. 
1. Physical Properties 


Refrigerant-oil nixtures have significantly 
different physical properties than pure refrigerants. 
Jensen and Jackman (Ref. 11) report that density and 
Specific heat behave ideally in refrigerant-oil mixtures, 
but that viscosity and surface tension do not. Ideal 


behavior of refrigerant-oil mixture density and specific 
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heat does not mean linear behavior. The governing equas 


tions are: 





1 2 E "€ 1 en C (2.5 
Pm Pol Pr2 
and 
E “05077: 0 
x ( ) GH 8 Sa (2-28 
where subscripts: 
- density ! = ከ ከ 
C = oil concentration mn = mixture 
Cy = Specific heat O- = TONI 
r = refrigerant 


Jensen and Jackman report that current refrigerant ص٦‎ 
mixture viscosity equations substantially underpredicted 
their experimental data: No predictive equation has been 
suggested for the surface tension of refrigerant-oil 
mixtures, thouyh Jensen and Jackman developed a correlation 
for EE 

Henrici and Hesse (Ref. 6] experimentally determined 
the surface tension for the R-114-oil mixtures that they 
used in their 1971 experiment. As shown in Fijure 2.329 BN 
surface tension of the nixture first decreased up to an oil 
concentration of 2.5 percent, and then increased continu 
ously with increasing oil Concentra oidis This type of non- 
linear behavior makes explaining the change in heat-transfer 
coefficient of refrigerant-oil mixtures, due to the changing 


physical properties of these mixtures, both difficulU 7002 
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possibly non-general. The behavior of refrigerant-oil 
mixtures is specific to the particular mixture components, 
and may be dependent upon the kind of oil being used. Some 
qualitative consequences of adding oil toa refrigerant, 
however, can be noted. 

The most observable result of adding oil to refrige- 
rants is foaming. Oil concentrations above 1 percent result 
in significant amounts of foaming from nucleate boiling. 
The foam bubbles form because the R-114ü in the R-114-oil 
mixture is more volatile than oil and vaporizes first, 
creating a gas bubble surrounded by an oil-rich layer (see 
Figure 2.4). Since the bubbles are coated with 011 2112 
with a higher surface tension than the bulk liquid anda 
have lower density, they rise to the top of the liquid. 
Because of their surface tension, the bubbles buiid up on 
the iiquid surface to produce a foam layer. 

This foaming action, which is most pronounced 
between 1-10 percent oil concentration (Ref. 6j, may affect 
the heat transfer of tube bundles significantly. For single 
tubes, it is the oil concentration gradient whicn would seem 
to play the major role, since the foam rises away from the 
tube surface and could only interact with the tube а 16 
sweeps by it from the bottom to the top of the tube. 

The general decrease in the heat-transfer coeffi- 
cient upon adding oil to pure refrigerants (recall Figure 
2.1) is subject to many explanations. Thome [Ref. 14], in 
an extensive review cf the literature, reports that the 
first explanation for the decrease in the heat-transfer 
Coefficient of mixtures was presented by Van Fijx et al. in 
19565 The effect was explained as being the result of the 
evaporation of the more-volatile components, leaving an 
oil-rich layer with a higher local boiling point, which 
increases the amount of superheat required to continue 


vaporization and bubble growth, thus reducing the 
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ecd 


dë 


763۰6-۲۶4257 ڈز6‎ ۲٢ ۶ ٣ Stephan and Preusser [Ref. 12] 
demonstrated conclusively that the work of formation of 
bubbles in a mixture is greater than in an equivalent pure 
2+2 They concluded that the mixture heat-transfer coef- 
ficient is lower than for tne equivalent pure fluid because 
of the resulting decrease in the bubble population. Àn y 
variation from a decreasing mixture heat-transfer coefficent 
(like Henrici and Hesse show in Figure 2.2) is attributed by 
Stephan [Ref. 13] to the non-linear variation of the phys- 
ical properties. Stephan proposes that the plot of surface 
tension (Figure 2.3), along with thermal-property varia- 
tions, accounts for the anomalous rise in the heat-transfer 
coefficient between 3-6 percent Oil concertration. 
Chongrongeong and Sauer [Ref. 10] state that it is the rate 
of heat diffusion, governed by the thermal properties of the 
oil-rich layer, that limits the bubble growth and that tne 
surface-tension effects are neglible. 

Thome (Ref. 14] proposes that all of the above 
factors, as well as the viscosity variation, are important 
in explaining the rise in the heat-transfer coefficient for 
sone refrigerant-oil mixture and heat-flux combinations. 

in summary, all researchers agree that the physical 
and thermal properties of a refrigerant-oil mixture are 
important factors in explaining the heat-transfer behavior 


cf mixtures. 


2. Saturation Temperature 


it has long been noted that increased saturation 
temperature ({i.e., increased boiling pressure) increases the 
boiling heat-transfer coefficient of surfaces in 
refrigerant-oii mixtures. In 1963, Stephan { Ref, 9) ٥ 
that at high oil concentrations, the heat-transfer coeffi- 
cient of refrigerant-oil mixtures becomes constant with 


respect to the saturation temperature. Stephan proposed 
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wide this because the addition of oil to a refrigerant 
introduces a large diffusion resistance, and that since the 
velocity of diffusion is almost independent of temperature, 
so should the heat-transfer coefficient become independent 
of temperature at high oil concentration. Figure 2:5 shows 
Henrici and Hesse's data on the effect of oil and saturation 
temperature in R-114. With no oil, the effect of raising 
the saturation temperature is seen to be a rise in the heat- 
transfer coefficient. With oil, raising the saturation 
temperature is seen to cause a slight drop in the heat- 
transfer coefficient. This effect has not been explained 


yet. 
3. Tube Dianeter 


Cornwell, Schuller, and Einarsson [Ref. 15] found 
that for smooth tube diameters from 6 mm to 30 nn, the 
nucleate pool boiling heat-transfer coefficient in pure 
refrigerants falls with increasing diameter. The effect of 


tube diameter was correlated by: 





Nu = ር peo (2:3) 
where 
Re - 4 D 
Һер y 
Nu = h D 
k 


Q 
ዘ 


150 for refrigerants 
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Since equation. (2.3) depends on physical parameters, the 
effect of oil could be significant. No known measurements 
iE effect of tube diameter Ino refrigerant-oil mixtures 


have, as yet, been performed. 


4. Surface Condition 


ua aa ee cmm: mto A A ي‎ EE vg 


Several researchers [Ref. 16, 17, 18] have investi- 
gated the e£fects of surface roughness on the heat transfer 
of pure refrigerants. As surface roughness increases, the 
heat-transfer coefficient vas found to increase due to 
increased nucleation. Nishikawa [Ref. 18] reported the 
effect of a variety of surface roughnesses in pure Е-114 
over a range of pressures from 0.294 MPa (42 psi) to 
2.94 MPa (820 psi). 


De Hydrostatic Effect 


el een A O AA ee = с ә 


The liquid cclumn above the boiling surface may 
generate large static pressures which will increase the 
oiling point. For R-114 at 0 OC (32 OF), a 0.3 m (1 ft) 
liquid head will raise the saturation temperature about 5 9C 
EXT). In large machines, this may be a significant 
effect. However, for small experimental apparatuses, the 


effect is negligible. 


B. NUCLEATE BOILING ОР REFRIGERANT-OLL MIXTURES FROM 
ENHANCED SURFACES 


Webb [Bef. 19), in an extensive review of the evolution 
of enhanced surface geometries, notes that the ability of 
roughness to improve nucleate boiling performance has been 
known for over 50 years. However, it vas not until 1968 
that the first commercial enchanced surface was patented 
[Ref. 20]. Since then, the number of commercial enhanced 


surfaces has dramatically increased as the understanding of 
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their design and operation has grown. Of the many possible 
methods for heat-transfer enhancement, tvo areas are 
currently being conmercially developed: 1) fins and surfaces 


with reentrant cavities, and 2) porous coatings. 
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Ihe ability of surface abrasion, open grooves, and 
fins to improve the heat-transfer performance of a smooth 
surface was first studied in the 1930's. The main dif mii 
culty with using surface abrasion to improve the heat- 
transfer performance is that fouling of the surface 
eventually returns the performance to that of a non-abraded 
surface. Studies in the 1950's and 1960's centered on fins. 
Recent comparative studies fcr refrigerants by  Carnavos 
(Ref. 4] and Yilmaz and Westwater [Ref. 2] found that fins 
and grooves result ina 50-100 percent permament improvement 
in the heat-transfer performance compared to a smooth plain 
tube in the same refrigerant. Webb [Ref. 19], in his liter- 
ature review, describes how researchers in the early 1960's 
found methods to improve the performance of fins by creating 
reentrant cavities on their surfaces. Reentrant cavities, 
such as shown in Figure 2.6, act as very stable nucleation 
sites and thereby enhance the heat-transfer perfcrmance. 
For a cavity to function as a nucleation site and remain 
active, even after the surface is subcooled, the mouth diam- 
eter (D) must fall within a critical range. Also, the 
Cavity must have a reentrant shape with a maximum reentrant 
angie (9). The optimum mouth diameter (D) and reentrant 
angle (6) are functicns of the fluid properties. 

The Gewa-T surface, patented in 1979 (manufactured 
by the Wieland Company), and the Thermoexcel-E surface, 
patented in 1980 (manufactured by the Hitachi Company), are 
two commercial surfaces which use modified fin shapes to 


form the necessary reentrant cavities. Figure 2.7 shows tne 
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Figure 2-6 Reentrant Cavity Geometry Factors. 
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(a) Schematic cross section of the Geomys s 
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Tunnel —” 


(b) Schematic view of the Thermoaex@e] = mom c 


Figure 227 Surface Details of Gewa-T and 
Thermoexcel-E Reentrant Surfaces. 
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details of these tube surfaces. Carnavos ( Ref. 4) found that 
in pure R-11, the Gewa-T surface outperformed a plain tube 
by 100-200 percent and the  Thermoexcel-E surface  outper- 
formed a plain tube by 300-400 percent. Gurvesüi0 o£ 
Figure 2.8 shov the relative improvement in the boiling 
heat-transfer performance of R-11 that can be achieved by 
using mechanically produced  reentrant cavities. 4ork 
continues to optimize these types of surfaces for the 
various refrigerants in use commercially. Both tubes have 
been tested in  refrigerant-oil mixtures and did not show a 
significant decrease in performance (Ref. 1 and 7]. The 
cost of these surfaces is not significantly higher than for 
smooth plain tubes, and the performance improvement is 


dramatic. 


2. Porous Coatings 


The second major type of enhanced surface is the 
porous boiling surface. Webb [Ref. 19] details the various 
production improvements and coating variations that have 
been made to the original 1968 patent by Milton of Union 
Carbide. The key to the performance of the porous coatings 
is their small reentrant cavities, which are interconnected 
by substrate tunnels. The particles used to make the coat- 
ings are usually copper or aluminum. According to Webb, 
researchers have found that the critical variable is the 
pore size rather than the particle size. Large pores are 
reguired for fluids with high surface tension aud high 
thermal conductivity. Small pores are optimum for fluids 
with low surface tension and low thermal conductivity (like 
refrigerants). Curve 6 of Figure 2.8 shows the relative 
performance of the High Flux surface to finned tubes and 
mechanically produced reentrant surfaces.  Carnavos [Bef. 4] 
found the High Flux surface to be 700-800 percent better 
than a smooth tube in R-11. No known studies have been 
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published on the performance of porous coatings in 
refrigerant-oil mixtures. Some studies of the neat-transfer 
performance of porous coatings in pure R-114 and 
refrigerant-oil mixtures have been made by Union Cartide, 


but their results are not found in the open literature. 
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III. DESCRIPTION OF EXPZRIMENTAL APPARAIUS 





A. OVERALL APPARATUS 


An overall schematic of the experimental apparatus is 
Shown in Figure 3.1, and a photograph is shown in Figure 
3.2.  Karasabun [Ref. 8] describes the design, Construction, 
and operation of the apparatus in detail. 

The apparatus consists of two Pyrex-glass tees. 2 ٦٣ 
R-114 is boiled in glass tee (1) and is condensed in glass 
tee (2). Gravity drains the condensate from the condenser 
back to the boiling section. A water-ethylene-glycol 
mixture at -17 °C (1 °F) 15 pumped through the condenser 
cooling coil via a computer-controlled valve (VC) to 
condense the R-114 vapor. The sump (7) that supplies tne 
water-ethylene-glycol mixture is cooled by a 1/2-Ton, R-12 
air-conditioning plant. 

Valve VC controls the R-114 liquid temperature and pres- 
sure. Figure 3.3 isa photograph of valve VC and the 
computer-controlled motor that operated VC. Opening VC 
causes more R-114 to condense and lowers the system pres- 
sure. Also, it returns more subcooled liquid to the boiling 
secticn which lowers the bulk liguid temperature. Singe 
data at many heat fluxes was desired for a constant tempera- 
ture, it can be seen that changing the heat flux without 
adjusting VC would change the system pressure and tempera- 
ture. A computer-controlled valve vas thought to be the 


best way to rapidly return the system to the desired satura- 


tion temperature following a heat flux change. Sections 
1350. 330011... describe. “im more detail the computer- 
contrclled valve and the operation f the system with the 


computer-controlled valve in use. 
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Photograph of Experimental Apparatus, 
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(b) Computer for Computer-Controlled Valve. 


Figure 3.3 Photographs of Computer-Controlled-Valve 
Equipment. 


oe 


Oil was added to the liquid 2-119 by draining it fro 
91. 2፡50) 00-7 š The oii cylinder was refilled as 
needed through valve V-2 from the oil reservoir (U). 

Two configurations of boiling tubes were tested. Short 
boiling tubes were tested to determine the correct assembly 
procedure to obtain data on the normal 431.8 un (17772 
long boiling tubes. The short tubes were cheaper to nake, 
thus more debugying attempts could be made by testing them. 
Section III.C describes the details of the construction of 


the short and long test tubes. 


B. CIL SAMPLING APPARATUS 


Following all data runs, an attempt was made to sample 
the local oil concentration in the vicinity of 8 bos ys 
tube. Figure 3.4 shows the oil sampling apparatus. BY 
opening valves S-1, 5-2, and S-4, the probe irne"camı በበኒ 
purged, trapping a sample inside a flexible s rilrcongas 
that was 30.5 mm (12 in.) long with a 3.16 um 1172717000 
inside diameter using pinch clamps. By weighing the sample 
tube and then boiling off the R-114 leaving behind the oil, 
the mass percent of oil in the R-114-oil mixture vas deter- 


nined. The mass fraction of oil was calculated by: 


Mass Fraction = mn. (320) 


m2 - ml 


where 
ml = mass of sanmple line 
m2 = mass of sample 11ре * Е-ТИФ икиси 
m3 = mass of sample line + oil (after boiling off R-114) 
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A Precisa Hodel 80 electronic mass balance was used to weigh 
the samples. The Precisa Model 80 is accurate: to 
7 0.0001 g. 


C. BOILING TUBE CONSTRUCTION 


1. Short Tubes 


Figure 3.5 shows the design of the short tubes. The 
short tubes were 15.9 mm (5/8 in.) in outer diameter, 
12.7 mm (1/2 in.) in inside diameter, and 203.2 mm (8 in.) 
in length. The short tubes extended 152.4 mu (6 in.) into 
the liquid R-115 from the left end flange. A 25.4 mm (1 in.) 
long epoxy plug insulated the right end of tne tube. እ 1 mm 
(0.04 in.) thick copper disk benind the epoxy piug was soft 
soldered in place to act aS a pressure barrier. The short 
tubes were heated by a 500-Watt 240-Volt stainless-steel 
cartridge heater. The heater was 6.35 mm (1/4 in.) in cuter 
diameter and 101.6 mm (4 in.) in length. 

The first short tube was made from thick-walled 
copper tubing. This tube was solid oxygen-free, high 
conductivity (OFHC) copper. Four 1.2 mm (3/64 in.) diameter 
holes were drilled into the wall of this tube at a diameter 
of 12.7 mm (1/2 in.) for thermocouple channels. Since this 
tube was solid, and had no sleeve interface, it did not have 
an interface resistance. Consequently, it was tne reference 
tube against which all other tubes were compared to deter- 
nine the amount of contact resistance they had. 

Six other short tubes were made. Five were mace of 
soft copper tubing and had sleeves inserted into the tute as 
indicated in Table 1. Soft soldering of the sleeves to the 
tubes was determined to yield negligible contact resistance 
by comparison with the solid tube. The last short tube (7) 
vas made of 90-10  ccpper-nickel and was coated with High 


Flux over the active 101.6 mm (4 in.) long section. This 
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(b) Thermocouple sleeve unwrapped (at section A-A) to show 
the relative locations of the thermocouple channels 
(all dimensions in millimeters). 





(c) Left-end view of the short tube. 


Figure 3.5 Sectional Views of Short Boiling Tube. 
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tube tested the compatibility of the soft-solder-assemnbly 
nethod with the copper-nickel tube and High Flux coating. 
Section V.A describes the results of data taken on the short 
tubes and the selection of the soft-solder method of 


assembiy for the long tubes. 


TABLE 1 
Summary of Short Tube Assembly Hethods 


Tube Surface Remarks 

1 Smooth Solid, thick-valled tube (reference) 
2 Smooth Slide-= It MON MIO ses 

3 Smooth Slide-fi1t (0.002 in. clearamec 

4 Smooth Press-fit (0.004 in. inter ose 9 

5 Smooth Press-fit (0.006 in. interference 

6 Smooth Soft-soldered 

7 High Fiux Soít-soldered 


2. long Tube 


Figure 3.6 shows the design of the long boiling 
tubes. These boiling tubes were 15.9 mn (5/8 in.) in outer 
diameter, 12.7 mm (1/2 in.) in inside diamter and 431.8 mn 
(17 in.) in llengtn: The center 203.2 nn (8 in.) vas the 
active test section. For the copper-nickel tube, the center 
secticn was the only portion of the tube that vas ccated 
with High Flux; The remaining 114.3 mm (4.5 in.) on either 
side cf the center section were smooth and unheated, and did 
not nucleate under any heat flux or oil condition. Karasabun 
(Ref. 8] describes how these  end-surfaces were treated by 
the data-reduction program as an extended fin from the 
center section and how their heat loss was accounted for. 

The center section was heated by a 1000-Yatt 
240-Volt stainless-steel  cartidge heater. The heater vas 
6.35 mn (1/4 in.) in outer diameter and 203.2 mm (8 in.) in 


length. The heater was surrounded by a copper sleeve with 


38 


'eqnj  6Purrrog 5uo7 yo MƏTA TPuUOT32ƏS 9”ዩ ,05ہ‎ 





1 
/ 
ኃበገሠቋ 11 ONIY-O //, 
2 


H31V 3H d Gang [ 1) — - —— 
3۸33 15 803033 


ና፣ 
እ ኣኣ Y REA DUNT 


.. 0 . 
E 244 Ким . 
Х Sasagot 6 a : 0 


111211712121121Ј111|11112111711111111111112711111| 111 7899950284899 





J3NNVH2 3811 ONI O8 
3100020۸111 | AE manne 


7 


7 
— ou 0-00 


39NV 14 ۸۷ 





/ 


// 


0 


eight 1.3 mn by 1.3 mm (0.050 in. 1777770 thermo- 


couple channels in them. Figure 3.7 shows the details of 
the channel layout. The thermocouple hot junctions were 
welded to the sleeve. Appendix A descrikes the calibration 


of the thermocouples. The channels were oriented to provide 
both axial and circumferential readings of the tube inner 
wall temperature. The sleeve was soft soldered to the tube. 
The data on the short soft-soldered tube closely matched the 
reference solid-tube data, and the short soft-soldered-tube 
data matched similiarly with long tube data. The maximum 
circumferential wall temperature variation in the long 
smooth tubes was 0.80 K (0.31 °F) at 50 k4/m2 compared to a 
solid tule circumferential variation of 0.34 K (0.61 LEE 
Section V.A describes in more detail the circumferential 
variation of temperature that resulted using the  varicus 
tube construction methods. Some axial temperature variation 
was experienced in the long tubes, particularly the long 
High Flux tube. Non-uniform heat generation from the 
cartridge heater is believed to be responsible for the axial 
temperature variation. Section V.B descibes in more detail 


the long tube axial temperature variation. 


D. DATA ACQUISITION AND REDUCTION 


A Hewlett-Packard 3497A automatic data acquistion/ 
control unit was used to read thermocouple outputs and to 
read an analog signal representing the current and voltage 
supplied to the cartridge heater. A Hewlett-Packard 98263 
computer unit vas used to control the HP-3497A and to 
analyze and store the data. 

Information was entered through the computer keyboard to 
prompt the HP-3497A to automatically scan each channel. All 
thermocouple measurements were accomplished by 0.245 mm 


diameter (30 gage) copper-constantan (type-T) thermocouples. 
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Bottom 


View of the boiling tube thermocouple locations as seen 
from the front of the experimental apparatus. 


Bottom 
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Thermocouple sleeve unwrapped (at section A-A) to show 
the relative locations of the thermocouple channels 
(all dimensions in millimeters). 
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End views of the boiling tube. 


Figure 3.7 Long Tube Thermocoupie Channels 
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A power sensing device, described by kKarasabun ۶٤٦ 
converted the ሏር current and voltage values supplied to the 
cartridge heater to a 0-5 V analog, DC signal for scanning 
by the HF-3197A. Table 2 lists the channel allocations in 
the HP-3497A. 

Following data acquisition for each point, results were 
computed according to the step-wise procedure outlined by 
Karasabun [Ref. 8], and summarized in Appendix B. Appendix 
B also includes a complete listing of the data-reduction 


program (DRP2). 


TABLE 2 
HP-3497A Channel Allocations 


Channel Purpose 

25732 Tube wall thermocouples ә. lon Tow 
Eo TIE) in ә tube, TI > J) nor 
used in short tubes) 


a 

Tí 
Liqui R-114 thermocouple T(9 
Backup liquid R-114 thermocou 
R-114 vapor thermocouple T(11 
Sump thermocouple T(12) 
Cartridge heater voltage analog signal 
Cartridge heater current analog signal 
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E. CCMPUTER-CONTROLLED VALVE 


The computer-controlled valve (VC) waS a Whitney, 
screwed bonnet, regulating valve with 3/8 in. Swagelok 
fittings. The valve travel was 10.5 turns from full shut to 
full open. The valve handle was replaced by a 101.6 om 
(4 in.) hard rubber disk (seen in Figure 3.3) which was 
rotated by the motor pinion gear. 

The computer-controlled motor vas a General Electric 
"Minigear Motor" with a speed of 105.7 rpm and a torque of 
3.39 N-m (30 Tb በከ The motor direction was controlled by 
two sets of Crydom solid-state relays that acted to open, 
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shut, or hold the motor depending on signals from the 
computer. The computer controlled the amount of time that 
the valve was moving. Approximately 3 seconds were required 
for the motor to turn the valve one turn. 

The flow through the valve was checked by a flow meter 
and verified to be approximately linear (1 turn - 10 percent 
flow) for the flows most often used. The valve position was 
tracked by a 10-turn potentioneter, which was connected to 
the computer-controlled motor by another 101.6 mm (4 in.) 
hard rubber disk driven off the motor pinion gear. Hard 
rubber disks were used instead of metal gears to avoid 
damage to the gear teeth during the program debugging stage. 
The rubber disks allowed sufficient slipping when, Lor 
example, the computer sent a valve-open signal even though 
the valve was fully cren. 

The R-114 liquid temperature input to the computer vas 
provided by a separate Copper-constantan  thermocougle 
installed in the same liquid R-114 thermocouple weli that 
the HP-3497A data acquistion/control unit used. The R-114 
theromocouple emf for the computer was amplified by an Omega 
thermocouple DC millivolt amplifier before input in to the 
computer. The data acquisition system andthe computer- 
controlled valve system were completely independent systens. 

The computer used for the computer-controlled valve vas 
an Octagon Systems SYS-2A microcomputer with an Esprit 1 
terminal connected via an RS-232C serial port. Appendix C 
lists the control program used. The control program was 
written in NSC "Tiny EASIC.' The control algorithm simulates 
a proportional-integral-derivative (PID) controller to vary 
the valve opening and shutting times. In the final program, 
limits were placed on the numerical value of certain program 
variables to prevent register overflow, jamming the valve 
fully open/shut, and to lessen the impact of system noise on 


the response of the systen. NSC Tiny BASIC is linited to 
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integers from -32000 to 32000 and fractional numbers OM 
truncated. The values of the weighting factors (A,B,C) for 
the prorortional term (E), integral term (I), and derivative 
term (D) were determined by trial-and-error. 


The control algorithm simplifies to the following lines: 


60 Input required temperature R 
230 Read R-1T4 temperature M 
260 Compute error n and change rin error (D) کے‎ 
262 Add 1 each loop to integral sum (T) lf error is positive 
264 Add -1 each loop to intégral sum (I) if error is negative 
400 Valve command V=(E/A) + (B*D)+ (1/C) 
420 If (V>0) then open valve 
430 Otherwise shut valve 
540 ССТО 230 

Section V.C desribes in detail the system response under 
this algorithm. Computer control of valve VC vas not used 
during the data taking as originally planned. Section 127 


discusses the reasons for this decision. 
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IV. EXPERIMENTAL PROCEDURES 


A. INSTALLATION OP TUBE IN APPARATUS 


Prior to installation in the apparatus, the boiling tube 
surface was cleaned with a 2 percent Nital solution (to 
remove surface oxidation and oil), rinsed with acetone, and 
air dryed. Short tubes tested with and without the above 
treatment showed no change in the heat-transfer perfcrmance. 
The treatment was effective in removing the slight surface 
oxidation present following soft soldering without changing 
the smooth or High Flux tube performance. Additionally, 
following testing with oil, a treatment was needed to return 
the High Flux surface to the "no-oil'" condition for further 
testing. 

After installing the tube in the glass tee, the appa- 
ratus was evacuated to 29 in. Hg by the portable mechanical 
vacuum pump (6) shown in Figure 3.1. System pressure was 
measured by a Marsh pressure gage (30 in. Hg to 150 psi 
mame, + 2.5 in. Hg and + 0.5 psi accuracy). The apparatus 


was left at vacuum for 2 hours to check for leaks pricr to 


each run. No noticable drop in vacuum was observed within 
the accuracy of the pressure gage. Next, the system gage 
pressure was raised to 0.19 MPa (27 psi), the saturation 


pressure of R-114 at 21 ОС (70 ОР), by opening valve V7 (see 
Figure 3.1 for the configuration of the valves). An 
Automatic Halogen Leak Detector, TIF 5000, was used to check 
ior R-114 leakage. The sensitivity of this detector is 3 
ppm minimum concentration. After pressure egualization with 
the R-114 reservoir (3), the reservoir drain valve (V6) and 
condenser return valve (V5) were opened to fill the boiling 
tee with liguid R-114. 
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Prior to installation, the left end flange had been 
Scribe marked to indicate the liquid level corresponding to 
1600 сс (2500 gm) of liquid RI4 aC 2I OI DE This 
vas the mass of pure R-11H48 in the apparatus at the beginning 
of the data” runs. The apparatus was now ready for taking 
data. 


B. GENERAL OPERATION 


Tabie 3 lists the 109 data runs accomplished during this 
thesis effort and their purpose. The data runs were 
numbered sequentially and preceeded by a ረ“ or 3-letter 


prefix to indicate the tube type. The tube prefixes were: 


Short Solid Tube 
High FIX ae Tube 
Wieland Hard 5 Long Smooth Tube 
Short Slide-Fit Tube 

Short Press-Fit Tube 

Short Soft-Soldered Tube 

Short High 10706 


The short tube runs consisted of 6 data points at 6 
different heat fluxes (usually 59, 37, 22, 14%, 070 
kh /m?). The normal tube runs consisted of at least 7 
different heat fluxes with 6 data readings at each heat flux 
(usually 98, 61, 37, 22, 14, 8, and 5 kW/m?) with some addi= 
tional low heat fluxes investigated for the O, 3, 7 
percent oil cases to check for the onset of nucleate boiling 
and hysteresis in the High Flux surface. 

In all cases, the data set was begun by starting the 
cooling pump (8) and opening valve VC slightly to slowly 
cool the liquid R-114 to the desired saturation temperature. 
The R-114 vapor temperature was also monitored and when both 
liquid and vapor temperatures were stable (usually after 
30-45 minutes), the heat flux would be established, the 
Saturation temperature reestablished, and the data aquisi- 
tion unit was allowed to take data. Usually, the vapor 
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TABLE 3 


Summary of Data Runs 


Remarks 
Runs WHO1 to WH10 taken by Karasabun 
Debug new tube f i 
PresS-fit tube 4 (0.004 in. interference) 


Effect of Tsat 
Tsat 
tube 

t 


2 (0.005 
Tsa 


in. clearance) 
Effect of | 
Effect of 7 


sat 
Rotate tube 90 degrees 


Study low heat flux error band 

E high heat flux error band . 
Rotate héater -90 degrees, tube fixed 
Shift heater and thermocoupies 
Repeatability 

Repeatability 

Repeatability 

Shift heater and thermocouples. 
Clean with Nital, acetone, and air dry 
0:0” 

Slide-fit tube 3 (0.002 in. clearance) 
Effect of Tsat 

Repeatabilit 


Tıme 55: Study (1 day later) 
Effect of Tsat 

Effect of Tsat 

Short Solid Tube 

557. 

Effect of Tsat 

4 days later 

Shift heater and thermocouples 

T days later 


Time Study Solid Tube (no effect) 
Short soft-soldered tube 6 
5257555” 

Debug short High Flux tube 7 

Study 757: superheat 

Repeatability 

Repeatability . 

Study hydostatic head (*1 in. level) 
-IncreaSe data pts at heat fluxes 
-Runs 49-54 form data set 

-Each run at different heat flux 
-Compare with set 55-60 

End of Set 49-54 

Study hydrostatic head 
-Same as above set 49- 
-same 

-same 


-same 

End ofwSet 55-60 ፣ 
Install Tvapor radiation shield 
Study Tliquid subcoolin 

Clean High Flux w/Nital and acetone 
Repeatability 

Repeatability 

Repeatability 

Rotate tube +90 and +180 degrees 
Repeatability 


ek in. level) 
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Table 3 


Summary of Data RUuS7C01u770) 
Run Isat No. 
No. (DCH Pts Remarks 
БИР 3 ee 7 “ez thermal grease to heater (no effect) 
SHF 70 6.7 36 Brfrect of 16501 
SHF71 -2.2 36 Re ھ2‎ i 
SHF72 -2-2 36 Debug ں ے 12ہ‎ 7٦ EIE CHE added). 
SHF73 -2.2 36 Remove and clean tube, add 1 percent Ga 
517۰7۷۲۲۰۳ 2702-6 2 percent oil 
591.) በመረ .ሪ NE 3 percent oil | 
Sür SDD SG 6 rercent oil 
SHr77 -2.2 36 “10 percent 0 . 5: 
WH78 E 9 Axial/Circumferential variation study 
ሽ፪73 = 2.2 BU? O percent oil, decreasinç q 
“H80 -2.2 11 Rotated tube -90 and -180 degrees 
WH8 1 -2.2 42 Repeatability GE n 7۶: 
WH8 2 ሮ ... ከ. percent omy ee i q 
14H83 -2.2 38 Repeatability OE Tun sS 
WH 84 6 ۷ 0 percent oil, increasing q 
WH85 22:2 2266 0 percent oil, decreasing q 
WH86 “22:16 1 percent oil, decreasing q 
WH87 6.7 48 1 percent oil, decreasing q 
WH88 2225055 2 percent oil, decreasing d 
WH89 6.7 48 2 percent oil, decreasing d 
VH90 2242080506 3 percent oil, increasing < 
WH9 1 22ے‎ MB percent oil, decreasing d 
WH92 6.7 ۵ 3 percent oil, decreasing q 
“493 6.7 38 6 percent oil, decreasing q 
WH9 4 ےب‎ 2۶ 6 percent oil, decreasing 4 
4Н95 -2.2 61 10 percent oil, increas mE] 
WH96 6.7 53 10 percent oil, increasing 3 
ИН97 -2.2 160 ٠٦0٣ ٣ 77 вәсл q 
VH98 -2.2 45 Repeatability (of run WH95) 
ዛ፲9 6.7 66 10 percent oil, decreasing q 
H510007 775 O Eercent oil, increasing q 
HETO SID 0 percent aill, 260 کہ‎ 
HF102 -2.2 24 Repeatability ርዕ. በ0051 27 0፻ 
፲፲ 1 0:3 6.7 66 0 percent oll, increasing q 
HF 104 6:7 639 0 percent. oil, E 
005 6.7 56 ۶ط 63۲4 ز6ط‎ ۱3 ٣ یں‎ ٦ u) 
HE GG ረ..ሻ 8 1 percent oll, decreasing q 
7000 ٦ 6.7 48 1 percent oil, decreasing q 
HF108 6.7 10 2 percent oil, decreasing q 
87109 መረ... 6, 2 percent oil, decreasing q 
HF110 6.7 48 3 percent oil, decreasing g 
НЕТТ 2762705” J percent oil, increasing q 
HEITI? ھچ ٔۀٴ‎ ٣ J percent oil, decreasing d 
HE113. —2.2 ٥٦ 6 percent oil, decreasing q 
HF114 6.6 48 6 percent oil, decreasing q 
HF115 -2.2 48 10 percent oll, nczes-175 74 
87116 -2.2 72 10 percent oil, decreas 
НЕ117 6.7 72 310 percent oil, decreas m m 
ba zk: 6.7 48 10 percent oil, increasing dq ` 
SPF119 -2.2 36 Complete study of time variation | 


Note: [ 
testing not listed. 
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75 Gata runs for Computer-Controlied Valve 


temperature would read higher, up to 2 K, than the liquid 
temperature due to the vapor becoming superheated. The 
apparatus, though cocled by the R-114, was still hot enough 
(60 OF) to superheat the vapor. Measurements of the liquid 
temperature, system pressure, and vapor temperature (with 
vapor probe shielded by a radiation shield) confirmed tnis. 
The liquid temperature best indicated the saturation 
temperature. 

The HP-3497A data aquisition unit voulü scan each 
channel, compute the heat-transfer coefficient, print the 
results (an example printout is shown in Appendix D), and 
store the data on the floppy disk. 

Folloving the taking of all data points, the data set 
vas statistically analyzed by subroutine STATS to compute 
the average heat-transfer coefficient at each heat flux and 
the standard deviation of the 6 data points for a given heat 
mU X. The standard deviation of the 6 data points was 
usually 0.5 percent for the heat flux and 1-2 percent for 
the heat-transfer coefficient. 

After all data sets at a given oil concentration were 
complete, oil was added via valve V1. The oil immediately 
dissolved in tne R-1184. No carryover to the condenser was 
noted except for several small drops at 10 percent oil ani 
the highest heat flux during the last few data sets. 
Foaming occurred with the addition of oil, and increased 
with both increasing heat flux and increasing oil ccncentra- 
tion. Figure 4.1 to 4.3 show photographs at heat fluxes of 
30 kWym? and 98 kW/m? for oil concentrations of 0, 3, and 10 
percent. When oil addition was not taking place, the oil 
cylinder and reservoir were isolated from the apparatus by 


valves V3 and V4 to minimize R-114 absorption by the oil. 
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(a) Oppercent oil ар kV/m”. No Foaming. 





(b) 1 percent Gilat 100 kW/m”. Foaming begins to appear. 


Figure 4.1 Photographs of Boiling and Foaming 
at 0 and Percent Oil. 
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s 3opercent oil at 37 kW/m*, 





(b) 3 percent oil at 100 kW/m2. 


Figure 4.2 Photoqraphs of Boilinq and F 1 
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(Ы) 10 Dercent Oil at 100 kW/m* . 


Figure 4.3 Photo papas of Boiling and Foaning 
at 10 Percent Oil. 
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C. SAMPLING OIL 


Oil sampling was done at room temperature to prevent the 


temperature of the sampling lines from boiling off the R-114 


and yielding false readings. The sampling procedure used 


vas: 
1. 


2. 


10. 


Weigh the empty flexible sample tubing and pinch 
clamps. Connect them to the sampling apparatus. 

Open S-1, S-4, and set S-2 to position 2 to provide a 
purge path for the sample flow. 

Lower the sample container to a height below the 
level of the R-114 glass tee. 

After purging the sample probe and lines, use finch 
clamps to isclate the 304.8 mm (12 in. flexible 
sanple tubing and trap a sample. Shut S-1 and switcn 
5-2 to _positicn 1. 

Disconnect the flexible sample tubing and immerse it 
in ice to lover its internal pressure below atmos- 
pheric pressure. 

Open S-3 to provide a vent for the purge container 
and pour the purge volume back into the boiling tee 
by lifting it to a height so that it flows by 
gravity. 

Weigh the sample line (the mass balance used was 
accurate to + 0.0007 g). 

Open a pinch clamp and allow air to warm the sample 
line and evaporate off the liguid R-114. Keer the 
Sample line on the mass balance to collect any drops 
of oil that may splatter as the liguid R-114 evapo- 
rates. 

After allowing several hours for the R-114 to evapo- 
rate, reweigh the flexible sample tubing. 

Calculate the mass percent of oil using equation 
13.41). 
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V. SYSTEM OPERATION AND PROBLEMS 


A. CIRCUMFERENTIAL WALL TEMPERATURE VARIATION 


Following the construction of the solid, thick-walled, 
reference tube, an investigation was begun to develop an 
assembly method which would result in negligible contact 
resistance between the inner tube wall and the copper 
Sleeve, and yield data comparable to that for the soiid 
tube. 

Short slide-fit tubes 2 and 3 (see Table 1) had sleeves 
that couid be slid into the test tube with tube 3 havinga 
tighter clearance. The slide-fit tubes exhibited a circum- 
ferential wall temperature variation of 7-14 XK at 3450 kW/m? 
in pure R-114 boiling at -2.2 °C. This variation matched 
data on a similar long tube (with slide-fit sleeve) tested 
by Karasabun [ Ref. 8]. Tube 2 had an average wall tempera- 
ture of 33 OC and tube 3 had an average wall temperature of 
25 °C at 50 kW/m2. The tighter clearance of tube 3 resulted 
ina lover contact resistance, and the wall temperature 
subsequently dropped. The short solid tube under similar 
conditions exhibited a 0.34 K circumferential wall tempera- 
ture variation and an average wall temperature of only 
19-7 7٦ Sauer et al. (Bef. 21), with a similar experi- 
mental apparatus, used a mechanically-press-fitted brass 
sleeve for obtaining data. Tubes 4 and 5 had copper sleeves 
that were mechanically cold pressed into the tube. The 
diametral interference was 0.01 mm (0.0004 in.) for tube 4 
and 0.015 mm (0.0006 in.) for БиБси The interface pres- 
sure obtained was calculated to be 15.2 MPa (2200 psi) for 
tube ü and 22.8 MPa (3300 psi) for tube 5. During rressing, 


the sleeves were lubricated with glycerin (C, H O, ) a 
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high-thermal-conductivity compound Елаз, accoradınq to 
Incropera and Dewitt (Ref. 22], should result in a 10-100 
times reduction in the contact resistance vhen combined vith 
the ccntact pressures mentioned above. 

Tubes 4 and 5 had circumferential wall tenperature vari- 
ations of 0.9-1.5 K with an average wall temperature of 
16.2 °C for tube 4 and 14.7 °C for tube 5 during initial 
testing. However, with tine, the circumferential wall 
temperature variation grew to 5-8 K and the average wall 
temperature increased to about 28 °C for both tubes. The 
resulting drop in heat-transfer performance is shown in 
Figure 5.1. This result is believed to be due to the phencm- 
enon referred to as "stress relaxation." 

Stress relaxation is a form of creep. The Metals 
decrease in stress resulting from transformation of elastic 
strain into plastic strain ina constrained solid. The 
phenomena occurs even at relatively low operating tempera- 
tures (at 25 °C, an 80 percent drop in stress can occur in 
200 hours), and is of most concern in applications like 
press-fits and solderless vrapped copper connectors. ihe 
time-dependent data of Figure 5.1 appears to be the result 
of a dropping press-fit interface pressure yielding a higher 
contact resistance, higher inner wall temperatures, and 
lower heat-transfer coefficients. 

Stephan and Mitrovic [Ref. 7] used a combination of 
mechanically press-fit and soft-soldered inner sleeves to 
obtain data on the Gewa-T surface in R-114. Soft-soldering 
was initially overlooked by this experimenter because the 
normal heating method used is an oxy-acetylene torch on the 
tube surface. The flame temperatures of an oxy-acetylene 
torch are above 800 9C (1500 OF). The resulting oxidation 
and heat damage (the High Flux coating melts at 800 9C) to 


the High Flux surface was unacceptable. Tests showed 
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however, that the inner cartridge heater could be used to 
apply a controlled amcunt of heat to melt a low temperature 
{melting point 160 °C) solder, while closely monitoring the 
tube temperature using the sleeve thermocouples. This 
method was used to produce soft-soldered short tube 6 which 
showed circumferential wall temperature variations of 0.8 K 
and an average wall temperature of 11.5 OC at 50 kW/m. 
Figure 5.2 compares the results of the soft-soldered and 
solid tubes. The agreement between the tubes is excellent. 
The long boiling tubes matched the short soft-soldered data 
very closely. 

Tests made by rotating the test tubes 90 degrees and 180 
degrees showed that the slight circumferential wall temera- 
ture variation of these tubes was due to the surface charac- 
teristics of the boiling tube rather than due to the 
thermocouples. Thermocouples located near more active 
nucleation sites of the smooth and High Flux surfaces had 
Slightly lower local wall temperatures. TROCOS OEK Circu- 
ferential wall temperature variation that resulted from 
soft-soldering the sleeves of the tubes is much smaller than 
the 11 K variation reported by Sauer et al. [Ref. 21], and 
is about the same as the 1 K variation reported by Stephan 
and Mitrovic (Ref. 7). 


B. AXIAL TEMPERATURE VARIATION OF LONG TUBES 


Karasabun [Ref. 8] reported an axial temperature varia- 
tion of about 20 K along the inner wall of his slide-fitted 
sleeve. Stephan and Mitrovic [Ref. 7] reported an axial 
temperature variation of 1 K along the sleeve they used. 
The long soft-soldered tubes tested in this experiment 
exhibited an axial temperature distribution that varied with 
heat flux. Figures 5.3 and 5.4 show the axial temperature 


ከ 5 ተሀህተ109 as a function of position along the boiling 
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surface for the smooth tube and Higk Flux tube respectively. 
The variation at 20 kW/m? (the heat flux at which Stephan 
and Mitrovic cite their 1 K variation) is less than 1 K for 
both tubes. The maximum axial temperature variaticn is 3 K 
for the High Flux surface at 98 kW/m2. 

The cartridge heaters used were precisicn-wound, 
сегапјс-соге, magnesium-oxide-insulated, Incoloy-sheathed 
WATLOW FIREROD heaters. The heaters had 6.35 mm (0.25 in.) 
long  lava-rock plugs at either end of the heater that 
reduced the actual heating length on each end. The heaters 
were initially believed to generate a uniform heat flux at 
all power settings, but the axial wall temperature data 
indicate the heat generation varied with the power level. 

Since the axial temperature distribution was fairly 
linear over most practical heat fluxes (less than 37 k1/n3), 
the arithmetic average of all 8 wall thermocouples was used 
to compute the heat-transfer coefficient. This results in 
Slightly lower heat-transfer coefficients and is a conserva- 
tive estimate of the performance of the High Flux surface. 
Appendix E analyzes the resulting error in the heat-transfer 
coefficient from using the arithmetic average for calcu- 


lating the heat-transfer coefficient. 


C.  COMPUTER-CONTROLLED VALVE 


Figure 5.5 shows the system response to closing valve VC 
manually and the improved response of the system when using 
the computer-controlled valve. By correctly cycling the 
valve between open and shut, the system saturation temrera- 
ture could rapidly be changed. Establishing a stable equi- 
librium temperature after a large valve movement was, 
however, difficult during the initial testing of the appa- 
Latus. The extensive trial-and-error testing for the 


correct weighting factors of the proportional, integral, and 
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derivative terms resulted in a large data base concerning 
the saturated system operating characteristics and a large 
amount of practice for the operator in manual control of 
valve ፻ር following 10 0 E UL computer-ccntrolled 
transients. 

Manual operation of valve VC was done in conjunction 
With the section of the data-reduction program (lines 1890 
to 2485 of program DRE2 listed in Appendix B) that monitored 
the saturation temperature and the rate of change of satura- 
tion temperature continuously. It was found that, because 
of the better temperature resolution of the HP-3497A data 
acguisition unit, with manual operation and a trained oper- 
ator, temperature transients could be maintained within 
OO This was a much tighter control bard than was 
possible with the ccmputer-controlled valve because of the 
poor resolution of its temperature sensing input from the 
Omega thermocouple DC millivolt amplifier. Additionally 
after much practice and experience, the manual method was 
found to be as fast as, or faster than the use of computer- 
controlled valve. 

To obtain the highest accuracy data inthe shortest 
period of time, computer-control of valve VC was abandoned. 
Section VII.B includes recommendations for further improve- 
ments to the computer control System to restore its 


usefulness. 
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VI. RESULTS AND DISCUSSION 


A. BOILING PERFORMANCE OF SMOOTH TUBE 


Figure 6.1 shows the nucleate pool-boiling performance 
or the smooth copper tube in R-114. The behavior wita no 


oil represents typical nucleate pool-boiling performance. 


Fcon point A toT polnt bB; a continuous increase in wall 
superheat (e D is observed when heat flax Is 
increased. No bubbles vere observed in this region of the 
curve as this region represents natural convection. At 


point B, incipient nucleate boiling occurs. From point eae 
point C (or C' for 10 percent oil), a reduction 117777 
superheat is observed while the heat flux is continucusly 
increased. This region is known as the aixed-boiling 
region, where transition from natural convection to nucleate 
pool boiling takes place. In this region, the heated 
portion of the tube began to activate an increasing number 
of nucleation sites, while the unheated ends showed no 
bubbles. In fact, the unheated ends underwent only natural 
convection, due to axial conduction of heat along the tube 
wall, at all heat fluxes. The transition from naturam 
convection to nucleate boiling occurred rapidly when there 
was no oil present. The surface would burst into nucleate 
boiling in less than a second after the first nucleation 
Site became active. At point C (or Ct), all ٤ ضط‎ 87 
nucleation sites were apparently active. After point C, the 
wall superheat again increases with increasing heat flux as 
Shown in region C-D. In region C-D, no new nucleation sites 
were seen to become active. Instead, the bubble departure 
rate increased. When the heat fiux is decreased after having 


established complete nucleate boiling, the curve follows a 
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different path (point D to polu. The existence of 
stable nucleation sites, which remain active over a wide 
Tange of heat fluxes, results in better heat-transfer 
performance than in natural convection, resulting in a lower 
wall superheat. 

The effect of adding oil is, according to Stephan 
[Ref. 9], to introduce a mass diffusion resistance and lower 
the heat-transfer coefficient. As seen in Figure 6.1, in 
region A-B, both the 3 and 10 percent curves were lower than 
the 0 percent oil curve. The 3 percent oil curve is lowes 
than the 10 percent oil curve probably because of the non- 
linear physical prorerty variations of Trefriqerant-ofri 
mixtures. The non-linear variation of surface tension (see 
Figure 2.3) would not seem to be responsible for this 
anomally. The curves in region A-B support the contentions 
of Chongrungeong and Sauer [Ref. 10] and Thome [&ef. 14] 
that the non-linear variation of physical properties of 
refrigerant-oil mixtures, other than surface tension, 
explains the heat-transfer behavior of Trefrigerant-oil 
mixtures. The effect of adding 10 percent oil was to delay 
the transition to ccrplete nucleate boiling on the tube. 
With 10 percent oil, the surface developed patches of 
nucleation sites that spread Slowly with increasing heat 
flux, until they covered the entire surface (point ከጠ 
Region Г-Е (or D'-E) shows that oil increased the wall 
superheat silghtly for 3 percent oil and significantly for 
10 percent oil. Again, this agrees with the concept of an 
increased mass diffusion resistance by tne addition of oil. 

Figure 6.2 shows the heat-transfer coefficient of the 


smooth tube in R-11lü-oil mixtures as a function of heat 


fl sx The curves show the heat-transfer coefficient of the 
smooth tube in region D-E, after complete nucleate boiling 
has been initiated. The effect of adding less than 6 


percent oil is seen to be small (about a 10 percent 
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reduction in the heat-transfer coefficient). However, an 
oil concentration of 10 percent causes a more significant 
drop in the heat-transfer coefficient (from 0 со 35 percent 
depending on the heat flux). 

Figure 6.3 shows more easily the degradation that cil 
causes in the boiling heat-transfer performance of the 
smooth tube. This figure plots the heat-transfer coeffi- 
Cient of the smooth tube relative to the heat-transfer coef- 
ficient in pure R-114 as a function of oil concentration. 


The effect of oil can be seen to depend also on the heat 


flux. Oil can be seen to generally degrade the performance 
of the smooth tube, except at heat fluxes less thar 5 kW /ne 
and oil concentrations between 2 and 8 percent. This 


behavior vas also seen by Henrici and Hesse (see Figure 
227: Since this curve shows the heat-transfer performance 
ın region D-E; with complete nucleate boiling, the non- 
linear variation of the physical properties of refrigerant- 
oil mixtures, including surface tension, again rorarii 
accounts for this ancmalous behavior. Since no measurements 
of the physical properties Of the R-114-o0il mixture vere 
made during this investigation, any possible non-linear 
property variations of the mixture used are unknown, 


requiring future work. 


B. BOILING PERFORMANCE OF HIGH FLUX SURFACE 


Figure 6.4 shows the nucleate pool-boiling performance 
of the High Flux surface in R-114-oil mixtures: The small 
magnitude of the wall superheats obtained during nucleate 
boiling should especially be noted. The High Flux surface 
in pure R-114 showed typical natural-convection region (A-3) 
behavior. The incipient point (B) occurred at much lover 
heat fluxes and superheats than for the smooth tube as shown 


earlier. The transition to nucleate boiling (B-C) was 
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sinilar to the smooth tube, occurring very rapidly (less 
than 1 second). Nucleate boilng fron the High Flux surface 
results in a drop of the wall superheat by about a factor of 
10. This is due to the extremely high density of nucleation 
sites present on the High Flux surface. From the low heat 
fluxes (1-10 kKW/m?) at which the transition to nucleate 
boiling occurs, it can be seen that the High Flux coating 
also assists in the activation of stable nucleation sites. 

Adding oil to the High Flux surface appears to delay the 
onset of nucleate boiling (point B wall superheat increases 
with increasing oil concentration). However, the transition 
to nucleate boiling still occurred very rapidly, even at 10 
percent oil. The rapid transition to nucleate boiling (less 
than 1 second) on the High Flux surface is probably due to 
the interconnected cavities which can assist in nucleating 
the entire surface once a single site becomes initially 
active. Oil is unlikely to inhibit this characteristic of 
the High Flux surface though it apparently delays the 
initial activiation of the first nucleation site. As seen 
in region D-E, for heat fluxes less than 37 kW/m2, the 
effect of adding oil to the High Flux surface, once it has 
been nucleating fully, is to cause about a 30 percent 
increase in the wall superheat. At heat fluxes in excess of 
37 kW/m* and at 10 percent oii, the wall superheat increased 
dramatically. As seen in Figure 6.1, the wall superneat at 
a heat flux of 98 kW//m? and 10 percent oil is about the same 
for both the smooth tube and High Flux surface. 

Figure 6.5 shows the heat-transfer coefficient of the 
High Flux surface as a function of heat flux for various oil 
concentrations. Again, oil is seen to degrade the nucleate 
pool-boiling heat-transfer performance. 

Air-conditioning plants typically operate in the heat 
flux range of 10 to 40 xW/m? with less than 1 percent oil. 


In this region of  pracical interest, the boiling 
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heat-transfer coefficient of the High Flux surface is about 
10 times better than a smooth tube. Also, it experiences 
only a 20 percent drop in the boiling heat-transfer coeffi- 
cient with the addition of oil (1 to 10 percent). This is 
seen clearly in Figure 6.6 which plots the heat-transfer 
coefficient, relative to the heat-transfer coefficient for 0 
percent oil, as a function of oil concentration. From 1 to 
10 percent oil at a heat flux of 14 kW/m?, the heat-transfer 
coefficient is about 80 percent of the no-oil heat-transfer 
coefficient. 

Figure 6.6 also shows that the oil-caused degradation of 
performance on the High Flux surface is nearly independent 
of oil concentration at practical heat fluxes. Only at a 
heat fluxof 98 kW/m? and 6-10 percent oil, does the High 
Flux surface experience significant performance degradation. 
At high heat fluxes and oil concentrations, the effect of 
oil may be to "clog" the interconnecting cavities of the 
High Flux surface due to boiling off of the R-114. Clogging 
the R-114 surface with oil would prevent replenishment of 
the nucleation sites with R-114 liquid, preventing the 
nucleation process and leading to higher superheats. The 
time-dependent behavior of the High Flux surface in high oil 
concentrations and at high heat fluxes was not studied in 


this experiment. 


C. COMPARISON OF BIGH  FLUX TO SMOOTH TUBE  BOILING 
PERFORMANCE 


rigure 6.7 shows the heat-transfer performance of both 
the High Flux surface and the smooth tube as a function of 
heat fiux. Again, the 7-10 times improvement in the heat- 
transfer coefficient by the High Flux surface is easily 
seen. At extremely high heat flux and high oil combinations, 
the High Flux surface performs oniy sligthly better than the 
smooth tube. 
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Figure 6.7 also shows that the heat-transfer curves of 
the High Flux surface and the smooth tube are not parallel. 
At low heat fluxes, the High Flux surface is more effective, 
in comparison to the smooth tube, in enhancing nucleation 
than at moderate heat fluxes. At high heat fluxes, the 
performance of the High Flux coating begins to converge 
toward the smooth-tube curve because the surface became 
vapor blanketed. This is expected since both tubes should 
perform about the same when they are completely vapor 
blanketed. 

Figure 6.8 shows the relative improvement of the High 
Flux surface over the smooth tube as a function of oil 
concentration. For the heat fluxes of practical interest in 
air-conditioning plants, the High Flux surface is 7-10 times 
better than the smooth tube. 


D. EFFECT OF SATURATION TEMPERATURE ON BOILING PERFORMANCE 


As reported by Stephan [Ref. 9], the effect of 
increasing the saturation temperature for both the smooth 
tube and High Flux surface was increased heat-transfer 
performance. Figure 6.9 shows the improvement in heat- 
transfer performance in pure R-114 achieved by raising the 
Saturation temperature? fron -2.2 1110171٦ to 6.19 ٣٦ 
(44 OF). At high heat fluxes, little improvement is seen in 
the High Flux surface performance because the surface is 
nearly vapor blanketed with bubbles. 

Figure 6.10 shows the effect of saturation temperature 
on R-115-oil mixtures. Increased saturation temperature is 
again seen to improve the performance of the High Flux 
surface as well as the smooth tube, even with 10 percent 
Oil. This is consistent with the no oil resis» but 


contradicts Henrici and Hesse"s data (see Figure 2.5). 
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E. LOCAL OIL SAMPLING EFFORTS 


The results of the attempt to sample oil locally near 
the boiling tube were inconcluSive. The method outlined in 
Section IV.C did not produce either repeatable or accurate 
results. Checks of the sampling method were made by 
sampling the bulk liquid with 0 and 10 percent oil at zero 
heat flux. The 0 percent oil check yielded oil concentra- 
tions from 0 to 2 percent. The 10 percent oil check yielded 
oil sample concentrations from 5 to 25 percent. 

Further refinements to the sampling apparatus are needed 
to permit an accurate local sample of oil in the vicinity o£ 
a boiling tube. Section VII.B contains recommendations to 


improve the sampling apparatus. 
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VII. CONCIUSIONS AND RECOMMENDATIONS 


CONCLUSIONS 


ከ. 


Are f Ill, the pool boiling heat-transfer coeffi- 
cient of the High Flux surface is about 10 times 
larger than that of a smooth tube. 

The High Flux surface began nucleate boiling at low 
heat fluxes, about 1 kW/n?, compared to a heat flux 
of about 10 k*5/m? for the smooth tube. 

Oil delayed the onset of nucleate boiling with the 
High Flux surface. 

Oil resulted in about a 20 percent reduction in the 
heat-transfer coefficient of the High Flux surface 
for most practical heat flux (less than 37 kW/m?) and 
oil combinations (less than 6 percent). 

At heat fiuxes of 98 kW/m? and greater than 6 percent 
Onl, the performance of the High Flux surface 
degraded by as much as 80 percent. The performance 
of tne High Flux surface was little better than the 
smooth tube at 98 kW/m? and 10 percent oil. 

The boiling heat-transfer coefficient of the High 
Flux surface was about 7 times better than that of a 
smooth tube for oil concentrations from 1 to 10 
percent over the range of heat fluxes employed in 
alr-conditioning plants (10-40 kW/nm?). 

The boiling heat-transfer coefficient of both the 
High Flux surface and the smooth tube increased with 
increasing saturation temperature. The improvement 
decreased at high heat fluxes for the High Flux 


surface due to vapor blanketing of the surface. 
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Local oil sample results were not accurate or repea- 
table enough to determine if a variation in the local 


oil concentration occurs near a boiling tube. 


RECOMMENDATIONS 


T 


The present cartridge heaters should be replaced with 
nore reliable heaters that will produce a uniform 
heat flux axially at all heat loads. 

A solid copper tube should be coated with High Flux 
and tested. Lue to the very smali superheats of the 
High Flux surface, even a small amount of contact 
resistance could affect the data at low heat fluxes. 
The physical properties of the R-114-oil mixtures 
tested should be measured to obtain information to 
better explain the reasons for the heat-transfer 
performance of both the High Flux surface and the 
smooth tube. Particularly, the anomalous rise in 
heat-transfer performance of smooth tubes when 1| to 6 
percent oil is present in R-114 should be studied. 

A secondary heater should be installed in the boiling 
section to keep the total heat input constant. The 
heat input to the secondary heater snould be 
increased when the neat input to the boiling tube is 
decreased and vice versa. This modification will 
maintain a constant heat load on the condenser and 
eliminate the need to operate valve VC, except to set 
the saturation temperature at the beginning of a run. 
Alternately, the effort to computer control valve VC 
Should continue by obtaining an accurate thernistor 
with a high temperature resolution for use with the 
SYS-2A mucroccHpulbems Improved temperature resolu- 
tion would allow the computer-controlled valve to 
operate properly and free the operator from the 


demands of manual control of valve VC. 
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The performance of a bundle of High Flux-coated tubes 
in R-114-oil mixtures should be studied. 

Data should be obtained on the High Flux surface over 
a wider range of temperatures and with oils of 
varying viscosities. | 

The High Flux surface should be tested for tinme- 
dependent  heat-transfer performance by taking data 
periodically over a long time during boiling. 

A small oil sample container with valves and a vent 
line should be manufactured. The flexible tubing and 
pinch clamps used in this experiment did not properly 
hold or vent the sample. 
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APPENDIX A 
THERMOCOUPLE CALIBRATION 


Karasabun [Ref. 8] describes the thermocouple calibra- 
tion equipment in detail. Two thermocouples were calibrated. 
One was made from wire at the beginning of the roll, the 
other from the end of the roll, following the making of all 
thermocouples used in the apparatus and tests. 

Essentially, the manufacturer-supplied calibration equa- 
tion for the thermocouple wire, a seventh order pclynorial, 
was corrected slightly by adding to it a second order curve 
fit of the variation of the manufacturer's computed tempera- 
ture for a given emf from a known set of reference tempera- 
tures (measured using a Hewlett-Packard 2804A quartz 
thermometer with a temperature resolution of + 0.0001 K and 
accuracy of trO OST]: 


The manufacturer's emf to temperature conversion egua- 


tion is: 

٣۰2 ۶۶ © .. * ase? * (እ. 1) 

a r” ተ a ES + a, ES + a, El 

where 

T = temperature (°C) 

ag = 0.100860910 

a, = 25727: 94309 

а, = =/ 07345-46295 

a4 = و7‎ 7٦1 

መ -92047436589 

ap = 6.97688E+11 

ac = -22606 3 

а, = 3.94078E+14 

E = thermocouple reading (volts) 
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Figure A.1 shows the quartz thermometer reading minus 
the thermocouple readings (discrepancy) versus temperature. 
The two thermocouples agreed to within 0.05 K of each other 
and the manufacturer's seventh-order polynomial needed about 


a 0.1 K increase to more accurately convert the enf's to the 


true temperature. The correcting second-order polynomial 
was: 
= + ወጪ T2 2 
DCP b. b T ot ( ) 
where 


DCP = discrepancy (K) 


እ. —292626899068E-2 
5. 2365099 07F- 3 
bə 52٣-5۲992698-5 
T = thermocouple reading 
(from equation A.1) (9C) 
Bus, the temperatures computed by the data-reduction 


program (DRP2) were emf's converted to temperature by equa- 
tion A.1 with corrections for that temperature computed by 
equation A.2 added to the temperature to get the true 
۲۰۱06 ۶۵3 ۲7۰. 

Since the data-reduction program utilized differences 
between thermocouples in all computations, such as wall 
temperature minus saturation temperature, the corrections 
above were necessary only for the computation of items 
dependent on the absclute temperature, like the fluid prop- 
erties. Appendix E describes in detail the uncertainty 
analysis and the effect of wall temperature variation on the 


computation of the heat-transfer coefficient. Thermocouple 
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variations of the soft-soldered tubes were most likely due 
to a slight amount of contact resistance and to the surface 
characteristics of the tubes tested, Since data runs 
involving shifting the wall thermocouples did not affect the 
data. 
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APPENDIX B 


DATA REDUCTION PROGRAM 


The data reduction program below consists of the 


following sections: 


Mai 
Sub 
Sub 
Sub 
Sub 
Sub 


n Program - 
Main - 
Flot - 
Poly - 
Plin - 
Stats = 


Menu of subprogram options 

Take data or reprocess data 

Plot data on log-log scale 

Ccmpute least-squares curve fit of data 
Plot data on eLinear-linear scale 


Ccnpute average and standard deviation oz data 


Subprogram Main consists of the following steps: 


1. 
SCH 
5 


Create data 
Select tube 


file for data and plot points. 
7 


Monitor heat flux or saturation temperature to estab- 


lish steady-state conditions. 


Scan all channels listed in Table 2 and save in data 


file. 


Ccnvert raw erf's to temperatures, current, and 


voltage. 


Conpute the 


hater. 


Compute tnhe 


heat-transfer rate for the cartridge 


average wall temperature, the wall super- 


heat, and the film temperature. 


Compute the physical properties of R-114 using given 


correlations at film temperature. 


Compute the natural-convection heat-transfer coeffi- 


cient of R-114 for the non-boiling ends of the test 


tube. 


10. Conpute the heat loss from the non-boiling ends. 
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11. 


Calculate the corrected heat flux 


to the liquid R-114. 


12. Calculate 


the boiling 


the R-114 fron the test tube. 


13. 


Print data. 


in plot file. 


The following is 


tion program 


Hewlett-Packard 9826 


a listing of the 


(DRP2) written in Basic 3.0 


computer. 
FILE NAME: 


DATE: 
PEV] SED: 


DREZ 
October 19, 1984 
March 102. 1985 


BEEF 

PRINTER IS 1 

FRING USING "4x ""Select ootion:^"" 

PRINT USING "EX “70 Jeking dala cer re-processing orevicus data' "^ 
PRINI USING ×58ٴ‎ "٠ Plotting data en Leg-Leo “"” 

FRIN! USING "6X,"72 FIetting data on Linear "^" 


PRINF USING *“6X,""2 Mate cross:olot coefft file””” 
۱٢۱۲۱۱۱ Ido 

IF Idor@ NIUN CALL Main 

IF ao THEN CALL Flat 

IE fdo=2 THEN CALL Plin 

IF Tdpe23 ¡HEN CALL Coef 

END 

SUB Moin 


COM ¿Gel Ce7y Ical 

DIM EnmfCI2) TCI2D Dlac6) U2at60,01a3€6) Doatb ) аб) Дак 6) ዞርህ6(8) 
DA!^ 0.10606091 25727.94369 -767*%45,8295 78025595.81 

пл Га -9242486549 6.92688869811,-2.Б6192Е915 3.940?7BE * 14 

REA Coed 
l'HINIFH 13 
(LEAN 709 


TAI 


BEEP 
INFUI CENTER MONI DATE AND TIME (MM:DO:1IH+MM:55)” Dates 
ONET 709; 7 ID" iate f 
ONIPIA 27034,7107 
ENTER 709:0ate$ 
PRINT 
PRIMI ° Menlh, idate and time :“ Dated 
FREN? 
PRINI USING *10X ,""NDIE: Lncogr am name : DRP27^77 
ВГЕР 
INFUI. "ENFER DISE. NUHPER?. pn 
PRINT USING "ISA." Dish number = '"" ያያ”)ዐዘቨ 
pi £ p 
INFUT CENTER. INPUT. ዘ('ዛቂ (ዐ፣* ንዐዕኳ4ጳ,!፣፡| 1[፪ )” Im 
RI[F 
INPUT "ENILR fHERHDCOUPLE. TYPE (O:NI V, 1*1 D) ` Ical 
IF Im=Q@ THEH 
BFLP 
Input "GIVE A WANE FOR THE RAW DATA FILE” 02 files 
PRINT USING “16%, ""New file name: "” [4ለ”)ሀ2 የ፤1[ሮፄ 
ዩ፻ ቦቦ 
11:۱۱:۱۱١۱ "Input SIZE Of | mE QUAL ,5 ۱ء‎ 
CRE ALE POAT O? filet امدرت‎ 
SSIGN ٥٥۱ 2۰نا‎ 1002 fitet 
ОРММт FILE tet ۱۱۶ ۱۱۴۰ FROWN 


DI filet- "DUMMY" 

CREATE RUAL DI files ++ el 
ASSICN @rtel TO MI Files 
QUIPUT 8 ilelsDate'b 
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from the test tube 
heat-transfer coefficient of 
Store heat flux and wall superheat values 


complete data reduc- 


for the 


1275 BEEP 
1252 INPUT "GİVE A NAME FOR THE FLOT FILE" P files 
1285 BEEP | 

1290 INPUT “INPUT SIZE OF FILE BDAT^,S1ze2 

1295 CREATE BD^T P files Su:e2 : 

1:00 ASSIGN @Plot TO P_files 


0209 6 
1310 INPUT "ENTER NUMBER OF DEFECTIVE TCS (QOsDEFAULT D" ,Idtc 


1315 IF Idtc-0 THEN 

1320 Udtcl=0@ 

1325 Ldtcd=@ 

1320 PRINT USING "16X,"^No defective TCs exist"^"** 


1335 “END” IF 
1342 IF Idtc=1 THEN 
1345, SBEEP 


1350 INPUT "ENTER OEFECTIUE TC LOCATION" ,Ldtcl 
[355 PRINT USING "16X,""TC 1s defective at location "" D'Late) 


1360 Ldtc2=0 


1265 ENO IF 
1370 IF IdtceZ THEN 
127 BEEP 


1380 INPUT "ENTER OEFECTIUE TC LOCATIONS" Ldtcl ,Ldte2 

1385 PRINT USING "18X,""TC are defective at locations "",0,4X ,0^ sLdtcl,Ldtc2 
1390 END IF 

1295 IF Idtc>2 THEN 


1400 BEEP 
1405 PRINTER 15 1 
1410 BEEP 


1415 PRINT “INVALIO ENTRY" 

1420 FPRINTER 15 701 

1425 GOTO 1305 

1420 ENO IF 

1435 OUTPUT @Filelitdtel Ldtcc 

1440! [m=] cotton 

1445 ELSE 

1450 BEEP 

14585 INPUT “GIVE THE NAME OF THE EXISTING DATA FILE” O02 files 
1460 PRINT USING ”16X,"”Dld file name: "",144":102. fileS$ 
1465 ASSIGN GFi1le2 TO 02 file$ 

1470 ENTER 8F11eZ2iNrun 

1475 ENTER 9File2:;0olg$ 


1400 BEEP 
1485 INPUT “GIVE A NAME FOR PLOT FILE” ,P files 
1490 BEEP 


1435 INPUT “INPUT SIZE OF FILE BOAT" 512ወ2 

1500 CREATE BDAT P file$ SizeZ 

1525 ASSIGN 8Plot TO P file$ 

1510 PRINT USING "165X,"^This data set taken on : ”"" 1]4ለ”!(0310ሀ5 

[515 ENTER $FileZiLodtel ,Ldaiel 

1528 IF Ldtel 0 OR Ldtcz^O0 THEN 

1525 FRINT USING *!8X.""Thecmocounles uere defective at [ocations: "" ,2( 30,453 0" 1 
Kare Lote: 


1550 END IF 
1535 ENTER 8۱۱۹۷۶۷ 
1540 END IF 


1545 PRINTER IS | 

15560 IF [m=0 THEN 

1555 DEEP 

1550 PRINT USING *4X ""Select tube tvoe””” 

1555 FRINT USING *6X,""0-*Smocth 4 inch Ref""** 

1570 PRINT USIHG "5xX,*"1*Smooth 4 1nch Cu (Fress/SI1de)* "^" 
1575 PRINT USING BX Ezo Solder 4 inch COS 

1580 PRINT USING “6X,°"3=Soft Solder 4 inch HIGII FLUX’ "® 
15ቦ5 PRINT USING "5X."*42Ujeland Hard B inch""" 

1530 FRINT USING ^8X,^"5:H[6GH fIUX B ٣ ٣٠۳ 

1595 PRINT USING “6X ,""6=GEWA-K 19 ۱۰ع‎ ۰۳ 

1600 INPUT Itt 

1605 IF Itt'h THEN 

1610 BEEP 

1615 PRINT "IHVALID ENTPY"* 

1620 (010 1560 

1625 END IF 

1620 OUTPUT 6Fileliltt 


BSS END IT 
1640 FRINTER 15 791 
1645 PRINT USING "IX "Tube Tvoe: Lt DO Nt 
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1650 
1655 
16501 
1665! 
1670 
1675 
1682 
1685! 
1590! 
legs 
1700 
1705 
1719! 
1151 
1720 
1725 
1720 
17351 
17401 
1715 
1750 
1755 
17601 
1765) 
1770 
Dus 
1780 
17851 
17301 
1795 
1800 
1805 
18101 
1815! 
1820 
1825 
1820 
1835 
1840 
1845 
1550 
1855 
1860 
1855 
1870 
1875 
1880 
1825 
1890 
1895 
1500 
1305 
1310 
bs 
1320 
1925 
1930! 
1955: 
1940 
1345 
] 456) 
1955 
1356 
1965 
1979 
5 


BEER 
INPUT “ENTER OUTPUT VERSION (@=LONG ,1=SHORT .2=NDNE >” Iov 


Di=Diameter at thermocounle positions 

DATA .0111125,.0111125,,0111125,.0129540,.012446,.0129540,.0100965 
READ Üləte) 

01=01а(114) 


D2=Diameter of test section to the base of fins 

DATA .015875,.015875,.015875,.0158241,.015875..015824,.01270 
READ D2a(*) 

D2-D2a(Itt) 


Diəlnside diameter of unenhanced ends 

DATAS 021.0127, 0121, .0132..9127..0132,.9111125 
READ Dial») 

Di=Dial ltt) 


Do=Dutside diameter of unenhanced ends 

DATA .015875..015875 ..015875 .015824 .015875..015824,.01279 
READ Doale) 

DosDoat Itt) 


L=Lenath of enhanced surface 

ን 1... .]015..101ር..1015..102015.:..03፡:ሙ,.፡ኢ02፡ሙ 0፡ፈዐጋሬ 
READ Lale) 

L=La(Itt) 


4 
Lu=Length of unenhanced surface at the ends 
1777775 705, 1051 02540 07525 0752 0752 
READ Lual») 
LusLua( Itt? 


KcusThermal Conductivitv of tube 
DATA 398 244 344,45 344.45 244 
READ Fcuata) 

Kcu=Kcua(1tt) 

A-PI*'Do 2-Di1^2)/4 

P-FI*Dc 

dəl 

S< =Q 

5 0 

5«5=0 

0 پ ہ5 
Repeat: !‏ 

IF Im=@ THEN 

ON FEY @,15 RECOVER 1870 

FRINTER IS ! 

PRINT USING *4X,"*SELECT DPTION""*" 
PRINT USING "86X,""OsTAKE DATA""" 
PRINT USIHG '6X,"^I«SET HEAT FLUX*"" 
PRINT USING "8X.”"2=SET Tsat ""” 
PRINT USING “4% ““NOTE: FEY @ = EÜSCAPE”"“ 
BEEP 

INPUT Ido 

IF Ido=0 THEN 2495 


LODF (0 SET HEAT FLUX 

IF Ido=1 THEM 

DUTPUT 709: "AR AF82 AL53 VRS” 

RF CP 

INPUT “ENTER DESIRED Ado” Dado 

PRINT USING "4X," "DESIRED Qda ACTUAL Qdo”"” 
Err=1000 

51:71:10 72 

OUTPUT 7091*AS SA” 


Ји 


Sum2@ 

FOR Ji=1 TO 5 

ENTER 703:E 

SumsSum*E 

NEXT Ji 

IF Isi THEN UoltzSum*5 

İF I*2 THEN ñmp=E 

NEXT I 
AgdosUolteAmp/OP1*02*L) 

IF ABS(Aado-Dada) Err THEN 
IF Aada>Dodp THEN 

BEEP 4080..2 

BEEP 4000,.2 

BEEP 4200,.2 

ELSE 

BEEP 25S0.. 
EEEP.- 250... 
ВЕЕР 250,. 
END IF 
PRINT USING *4X,MZ.230E ,2X MZ.30E " iDqdo .Aqdo 
WAIT 2 

60T0 1970 

ELSE 

BEER 

PRIMI USING '"43X.H7.3DE , 2X, M7. DE " «Dado Aado 
0ع‎ ; 

VAIT 2 

60፲0 1975 

END IF 

EHO IF 


КӘ r r3 


LOOP TO SET Tsat 
16108221 


BEEP 

INPUT “ENTER DESIRED ٢ , 8۹ 

PRINT USING "4x,"^ DTsat ATsat Rate Tv 
Old1:0 

0142=0 


OUTPUT 7091۰۵8 ۲5ح۸ 37ھ‎ ٦٠ 
FOR I=I TO 3 

Sum-Q 

DUTPUT 708: “AS SA" 

FOR Ji=l TO 2 

ENTER 703:Eliq 

Sum=Sum+tElía 

NEXT li 

Elio=Sum/2 

TIid=FNTvsviFlio) 

IF Tel THEM AtldeTld 

JF 1=3 THEN 083ء۲‎ 

MEXT 1 

IF ABS(AtId-Dt!d>)>.2 THEN 

IF Atid Otild THEN 

REEP 4000,.2 

BEEP 4000..2 

РЕЕР 4202,.2 

ELSE 

BEEP 250.. 
DEER 250,. 
RECEP 259 
CHO If 
Erri=Atid-Old!l 
Oldl=At1d 
Err“1.-0Uldz 
Dld2"Ty 


fra r rə 


PRINT USTNG "4x,5(40D.DO 1x?*:0t 1d Atld Erri.fy Err? 
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2279 
2280 
2385 
2530 
2335 
2420 
2425 
2410 
2415 
2420 
2425 
2420 
2435 
2440 
2445 
2450 
2455 
24501 
2465 
2470 
2475 
2482 
23851 
2430! 
2496 
2500 
0> 
2516 
2515 
2520 
2525 
2520 
2535 
2540 
2545 
2550 
2595 
2558 
ፎ585 
2570 
2575 
25680 
2585 
2590 
2595 
2622 
2505 
2610 
5655 
2670 
20.3] 
26301 
2635 


WAIT 2 

6010 2170 

ESSE 

IF. ABS(AtId-Otld}>.1 THEN 
IF Atid»Otid THEN 

ВЕЕР 2000,. 
BEEP 3000 ,. 
ELSE 

REEP 800,.2 
ВЕЕР 820,.2 
END IF 
Егг1=А(14-0141 

Didi-Atid 

Err2=Tv-D1d2 

Old2=Tv 

PRINT USING "4X,S(HD0.DD ,1X)”iDtl1ld Atld Ercl Tv ,Err2 
WAIT 2 

GDTO 2179 

ELSE 

BEEP 

ErrisAtid-Old!i 

Üldlə3tld 

Err2=Tv-D1d2 


rJ ^J 


0122=Гу 


ሺ] 1 USING AX SI MDO LO |[ደ)7”1311/ AIG ECE TV Erz 
WAIT 2 i 
GOTO 2179 

EHO IF 

END IF 

END IF 

ERROR TRAP FOR Ido OUT AF BDIINOS 

İF 1d3o“2 THEN 

REEP 

GOTO 1890 

END IF 


TAKE DATA IF Im=@ LOOP 

EFEP 

INPUT “ENTER BULK DİL 7“ Boo 
OUTPUT 7069۱۰۵۸۳ ۶29ھ‎ AL35 URS" 
FAR I=! TO 12 

QUTPUT 7991"AS SA” 

۸۵9 ؟ 

FOR ,11*] ፲0 2 

ENTER 709:£ 

SumsSum*E 

NEXT Ja 

Enf CI )2Sum/2 

NEXT 1 

DUTPUT 709 “AR AF62 ALB3 URS” 
FOR I=! TO û 

QUTPUT 709:"AS SA" 

Sum=Q 

FOR Jı=1 TO 2 

ENTER 709:E 

SumsSum*E 

HEXT Ji 

IF (el THEM Ur-Sum/2 

IF 1:2 "HEN Irc*5Sum/Z 

HEXT | 

ELSE 

ENTER Fi l1e2:8oo .Tolds Emft(ə ) Ur Ir 
END IF 


COMERT emf'S TO TEHP ,UOLT CURRENT 
Tua= 0 
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2640 
2845 
2650 
2655 
2660 
2665 
2670 
2675 
2680 
2685 
2690 
2695 
2700 
2705 
2710 
2716 
2720 


57^ 
Lía 


2730 
205 
2740 
2745 
2750 
2455 
760 
2265 
2770 
2775 
2780 
creo 
2740 
2795 
2800 
2825 
2810 
2815 
820 
2825 
2830 
28255 
2840 
2845 
“850 
255 
2860 
2865 
2870 
.875 
2880 
2895 
28901 
2895 
2900 
305 
2910 
5915 
2920 
2925 
2920 
። 5 
23421 
2 
2450 
2 5 
2360 
2965 


FOR I=1 TO 12 

IF Idtc;0 THEN 

IF I-Ldtc! ОН 1422 ТНЕМ 
[(] )።>99, 99 

GDTO 2710 

END IF 

END IF 

IF 1tt-4 THEN 

IF 1>4 AND I<9 THEN 
Hi” 

GOTO 2712 

END IF 

END IF 
T(1)=FHTvsv(Emf(I)) 
NEXT 1 

IF 1tt<4 THEN 

FDR I=1 TD 4 

IF 1eLdtci DR I=sLdtc2 THEN 
Tua=Twa 

EL SE 

Twa=Twa+T(I) 

END IF 

MEXT 1 
Tua=Twa/(4-ldtc) 
ELSE 

FOR I-1 TO 8 

IF IsLdtc! OR I-Ldtc2 THEN 
Tua=Tua 

E1 SE 

Twa=fwatl( 1) 

ENO IF 

NEXT I 
TuasTua/(8-Idtc? 
END IF 

T1d-T(3) 

TlazəeT(10) 

TvsT(11>) 

IF Itt 3 THEN 

11 == 

Tvət T( 103:TC(11))/2 
END IF 

Tsimmo=T(12) 

Amo=1r 

۷۰ ٤ 

۱0-۷ Amo 

IF Itt=Q@ МЕН 
VcusFNEcutTwua?) 

ELSE 

F"cu*Fcuat( Ttt? 

END IF 


FOURTER CONDUCTION EQUATION WITH CONTACT RESISTANCE NEGLECTED 
TweTwa-OeLOG(D2 ‘DI IC fF i ekcuel) 

Thetab=Tw-T1+3 

IF fhetab @ THEN 

BEEP 

INPUT “TWALL TSAT (@=CONTINIE, 1-ENÜ)” Te, 

IF lev=Q4 THEN GOTO 1875 

IF lev=1 THEN 3330 

ENO IF 


COMPUTE "ARIDI!'S PROFFRTIES 
Ifilm=fHIfilm( Tu T Id) 
۲٦۷۶۲۲۱۴۶۰۰۰۲) ۲٢۱ ۸ا‎ ( 

MueFfiMu( Tfilm) 

۲٢٢۲٢۶٢۳٢۰۱۲۶ ۱ ۸ ( 
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2970 CosFNCotTfi1lm? 

2375 Beta=FHBeta(Tfilm) 

2980 Ni=Mu/Rho 

2995 AlohazK/(RhosCo) 

2990 PreNi/Alpha 

2995 Psat=FNPsatíTld) 

30601 

30051 COMHFUTE NATURAL -CÜNVECTUE HEAT-TRANSFER CDEFFICIENT 

3010! FOR UNENHANCEN END(S) 

3015 Hbar=190 

2020 Fesz(Hbar*P/(Kcu*A))^.5*Lu 

3025 Tanh-FNTanh(Fe) 

3020 Theta=Thetab*lanh/Fe 

3035  X.-(9.81*«Beta*Thetabs*Oo 3*«Tanh/(tFesNi*Alpha2) .166567 

3040 Yya(14(.553/Pr17(8/1581) (8/27) 

2035 Hbarc=kK/Do*(.B+.387*x+./Yy)"2 

3050 IF ABS((Hbar-Hbarc)/Hbarc)^.001] TIIEN 

3055 Hbar=(Hbar+Hbarc)».5 

3060 (6010 3070 

2055 END IF 

10701 

7075! COMPUTE HEAT LOSS RATE THROUGH UNENHANCED END(S) 

iQga Qls(Hbar*P*Ecu*A) .S«Thetab*Tanh 

2005 Qc40-2*01 

2090  As-PI*D2*L 

3095! 

2100! COMPUTE ACTUAL HEA) FLUX AND 8OILING COEFFICIENT 

3105  QdosQc/As 

2110 Htube-fğdo/Thetəb 

31151 

3120! RECORD TIME OF DATA TAKING 

1125 IF Imsà THEN 

3130 ООТРУГ 709:*10" 

3175 ENTER 7?Q091Told$ 

3140 END IF 

3145! 

3150! OUTPUT DATA TO PRINTER 

3155 FRINTER IS 701 

3160 (1 JovsQ THEN 

2155 PRINT 

3170 PRINT USING "“10X,* "Data Set Humber = “° DDD -X.""Bul* 011 % = "" DD.D,5X,I 

ЈА"*(Ј,Воо,12149 

3175 PRINT 

3180 PRINT USING ^10x,**TC No: | = 3 4 E 6 d 
gre. 

SES BRIT USING SiOx “Ten” BCIX MOD. OO TO), T(29.T631,16(4) T(5),T(89.T( 

0۰۰۱ 

3190 PRINT USING "1QX.,"^ Tura lliod Tlıqd2 T.əor Psat Tsunp” ”” 

2125) USİNE “19X 2CMDD.DO. 1x2). X HDD.OD 1X ,2C€CTX ,HDD.DD),2X , HDD.O" sTua,TId, 

Tid2,Tv,Psat,Tsump 

3200 PRINT USING “10X.7” Thetab Htube Qde” ”” 

3205 РР1МТ USING "10X ,MDD.30,1X .M2.3D0E, 1X , M2. 3DE" Thetab .Htube ,Qdo 

3210 END IF 

3218 IF lov=! THEN 

Vn a IF J=1 THEN 

ә PRINT 

32306 PRINT USING "10X."^RUN No 01l% Tsat Htube Qdo Thetab””“ 

3225 END IF 

3240 PRINT USING “3¿X,30 4x DD,2X.M00,D0,3(1X,M2.3D£)”"1J,80p .Tld Htube do ,Thet 

ab 

3245 END IF 

31280 IF Im=0 [HEN 

3255 BEEP 

3250 1ዞህ[ "OK TO STORE THIS DATA SET (t=¥ Q=N)7" OF 

3265 END IF 

3270 IF Dlel OR Ims! IHEH J51*1 
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275 IF Ok=l| AND 1:0 THEN QUTPUT 8Fi11el:Boo.ToeId$ .Emf(C*),Ur,ÍIr 
3280 IF Im=1 OR O=I THEN QUIPUT 6PlotrOdo.Thetab 
3285 IF Im=@ THEN 
2230  &EEP 
3295 INFUT “WILL THERE BE ANOTHER RUN (12Y,0-N)?" Go on 
3200 Nrun=J 
3305 IF Go on«^*| THEN 32370 
3310 IF Go_on=! THEN Reoeat 
3318 ELSE 
3320 IF J:Nrun+1 THEN Reo=at 
3335 ٣ 
3230 ١۱۴ 1:0 THEN 


3335 BLEP 
3240 FRIMT USING "1Qx ,""HOTF: "",22,"* data runs were stored in file ^^,106571J- 
٠۱۰۴2۶٤ = 


3343 ASSIGN Gf ilel TO * 

3750 QUTFUT 8File2iNrun-l 

3355 ASSIGN 9Frlel TO OI filet 

3250 ENTER @FiletiNatet® Ldtcl tdtcc Itt 
3365 OUIFUT ƏFilezZiDate$ Ldtci Ldte2,1tt 
3270 FOR I=1 TO Nrun-1 

2278 ENTER éFileliBoo .To!ldS Emf( +) Ur Ir 
37380 QUTPUT @File2:Bop ,TolJdé Emf( +) Ur,lr 
3385 NEXT I 

3390 ጸፍፍ16ከ 81121] [0 * 

33895 PURGE "OUMMY" 

3400 EMO IF 

2405 BEEP 

1310 PRINT USING "10X,""NOTE: “5.22.” X-Y oaucs were stored in plot data file 
271900 (0-1 Ele 

1115 ASSIGN @F:le? TO ə 

3420 ASSIGN @Flot JO + 

3325 CALL Stats 

3430 BFEP 

ux E INPUT SIE RESTO PLOT DATA: (TY: Q= Ni OF 
3442 ፲፻ 0፥።=1 THEN 

3445 CALL Plat 

3450 END IF 

3455  SURFNI) 

3460! 

324651 CURVE FIIS QF PROPERTY FUNCTIONS 
3470 DEF ЕНКси(1) 

3475! OFHC COFPER 250 TO 308 vr 

8+ + S IC TÜ 

3485 ኣ።=4:/4-.1125!ነ 

3490 RETURN K 

2435 FNEND 

3600 OFF FNMuCT) 

3505! 170 TO 360 K CURVE FIT OF UISCOUSII* 
3510 Tk=]+273.15 Li TO E 

3818 Mu*EXPC-41.46528*C1011. 37/7 TIE 2))* 1, OE-3 
2620 RETURN Mu 

23625 FNENO 

+0 0 7٦ 

36301 180 10 400 F CURVE FIT OF Co 

2540 TksT1273.15 Ic Tn k 

3545 CÜpə.40188:1 .65007E -3əfrl 1.51433E-BeTi “2-6.57952£-10əTk 3 
2550 CozÜos 1060 

3555 RETURN የፀ 

2262 FMEND 

ፓ565 DEF ፻| ቦክከዕ(!) 

32970 Ti =[+272.15 iC TA + 

5575 ጸ#፡1-(1,851) /7ጄ2,. 955) İK TOR 

3500 Ro=26.32+61.146314+ (| 2)፥16,4160በ15፡፪፥1!1?7,ፈ476626»፪”",5+|. || 9828»2”2 
"585 Ro=Ro*.0624723 

2590 RETURN Ro 
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3995 
3620 
3605 
3612 
3615 
3620 
BES I 
3630 
3635 
3540 
3645 
3650 
3655 
3660 
2665 
2670 
abb 
JEFO 
2685 
2590 
2695 
3700 
3705 
3710 
3715 
3720 
3725 
3720 
3735 
2746 
3745 
3759 
3055 
3750 
3765 
3770 
205 
3780 
3725 
3790 
3145) 
2600 
3806 
2810 
ə 
3020 
3825 
2810 
2029 
3840 
2045 
2960 
3855 
3РБО 
3865 
3870 
777 
2090 
1885 
3840 
3895 
3300 
3905 
3910 
Fan 
3920 


ЕМЕМО 

DEF FHPrCGT) 
Pr:FMCpCOCTO)«FNMuCTO/FNK(COT) 
RETURN Pr 

FNEND 

DEF FMK(T) 

TS380 K WITH T TN C 
K=.,071-,000261+T. 
RETURN K 

FNEND 

DEF FNTamhiX) 
P-EXP(X) 

0=1/Р 
Tenh-(P-Q)/(P30) 
RETURN Tanh 

FNEMD 

DEF FNT+s.(U) 

COM real 


۲:۱۱۵ ( 
FOR 10 7 
T-T4CCT oU" I ]آ‎ 
HEXT I 


IF Ical=1 THEN 
T=T-6.,7122¿934E-2+1+(9.0277043E-3-T+«(--9.3:53917E-5)) 
ELSE i 
T4T*8.625897E-2*T*( 3.76199E - 3- T«5.0583259E-5) 
END IF 

RETURN T 

FMEND 

DEF FNBetatT) 

ዞርፀ=ቶከቦክዕ(ፐ+፥. |) 

RemseEHPho(1-.13 
Beta=-2/(RootRomde(Rov-Rom)/.2 

RETURN Beta 

FMENO 

DEF ۲٢۲۲۶ 1۸۸٢) ٢سب‎ Ild) 

TfilmstTusTId'7/2 

RETURN Tfilm 

FNEMD 

CFF FHPsattTc) 

TO 80 deg F CURUE FIT OF Pset 

Tf*1.8*Tc*32 
Pa*5,9415525*Tf*(,153520R2 t [f£ *(1.484091653E-34 Tf *9.8150571E- 862) 
Вр=Ез-14.7 

IF Pq Q THEN I +=PSI6 ,-=in lig 

Psat=Pa 

ELSE 

Psat=P9:29.9:“14,7 

END IF 

PETURN Psat 

F NENN 

SUR Plot 

COM /Colv/ ACI@,.10) C(10) BC4) Nop ,Iprnt ,Opo [log 
INTEGER Ta 

PRINTER TS 1 

BF[ P 

INPUT “LIKE OFFAUI T UALUF?2 FOR PLDT (leY ,0:N)?” Tdv 
ВЕРР 

FRINT USING "4x "Galart Mption:””” 

РАТМТ USING “6X *"@ о versus delta-T""" 
PRINT LISING "5X,"*] h versus delta-T1^""* 
FRIMT USING "BX,""72 bh versus n'^" 

INPUT Ana 

ВРЕР 

IHFUI "SELFCT UNITS (@=S[ l1=FHGLISH)"” [un 
PRINTER TS 705 


ay 


2425 
3930 
2925 
39456 
2345 
3950 
35955 
2060 
3955 
2370 
3975 
39R0 
3985 
399A 
3935 
3000 
4008 
4010 
4015 
4020 
4025 
4070 
4235 
4242 
4045 
4060 
4055 
4060 
1025 
4076 
4075 
40RA 
4025 
4282 
4005 
4100 
4105 
4110 
4115 
3120 
4175 
4170 
4135 
4140 
41465 
4150 
4155 
4162 
4165 
4170 
4175 
4180 
4125 
4190 
4195 
3700 
4205 
4210 
4215 
$270 
1225 
1970 
4238 
4240 
4245 
4250 


IF 1dv:-1 THEN 

BEEP 

IHFUT "ENTER NUMBER OF CYCLES FOR Y-AXIS” Сд 
BFEP 

INPUT “ENTER NUMBER OF CYCLES FOR Y-AXIS" Cy 
PEEP 

INPUT "ENTER MIN X-VALUE (MULTIPLE OF 109)” Xmin 
BEEP 

INPUT “ENTER MIN Y-UALUE (MULTIPLE OF 109)” Ymin 
EE 

IF Ooo=@ THEN 

Cy=3 

Cx 23 

Xminz.]l 

Yminsz100 

EHD IF 

IF Ooo=1 THEN 

Cy=3 

C«23 

Xminz.] 

Yminz100 

END [F 

IF Ope=2 THEN 

Cv=3 

Се=2 

¥min=] AAA 

Yminz 100 

ENO IF 

END (F 

ВЕРР 

PRINT *164:5Р141Р 2322 2202 9200 .6800:” 

PRİNT "ፍር 4,109 A ICA1TL 2 M:” 

Sf.=1890/Ca 

Cfv=100/Cy 

PRINT "PU oo Pü: 

Nne8 

FOR {=f TON Court 

Yarzımınə 10”(1-1) 

IF [=ር.+1] THEN Nn=1 

FOR (ef TÜ fin 

IF J=1 THEN PRINT "M 2 0° 

1F J=? THEN PRINT “TL 1] በ” 

Xa-£ate] 

Yer GT(Xa/Ymın )əSf, 
PRİNT "РА":Х," 24 XTi” 
NEXTJ 

NEXT | 

PRINT “PA 100 APUI” 
PRINT “PU PA ñ ..0 PD" 
ኮ1በስጎ ። 3 

FOR I=1 TO (y+1 
Yat=Ymin+ 10" (T-)) 
IF I5Cy*«1 ፐ| ይክ ከክ።! 
FOR Jz) ፐበ ከክ 

IF da THEM PRINT "TE 2 0” 
IF 1-2 TIHEM FRINT “TI | 
va=Yato] 

Ys| GTC Ya/Ymin sS fy 
FRİNT “PA A “Y "YT" 
NEXT J 

HEXT 1 

PRINT "PA AINA TL 0 2“ 
Нл-8 

FOR Tat TAO ር». ፥1 
Xat=Xmine1iH (1-1) 

IF I26«41 THFN Мәз| 
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4255 FOR J=1 TO Nn 

4260 TIF (el TIIFN PRİNT “TL @ 2" 
4288 IF J`*1 1 THEN PRINT “TL à I" 
427 Xa=XateJ 

4275 X=1 GT(Xa/Amin)eSfx 

42889 FRINT "PAT," 10A: XT" 
425 NEAT 1 


2401 NEXT J 
4295 FRINT "PA 100 160 PU PA 100 9 PD” 
4300 ٦ ۹ 


285 "08 1መ=1 IO Cv*l 
4316 ሃፅፊ!=ሃጠ14በ*#|0” (1-1) 
4315 IF I=Cy+1 THEN Nn=1 
4320 FOR J=1 TO Nn 
4325 IF Jet THEN PRINT "TL @ 2° 
4320 IF J>1 THEN FRINT “TL I“ 
4335 YazYated 
4340 Y=lG6T(Ya/Y¥mindeSfy 
4345 PRINT "PD FA 100,^ Y "YT” 
4350 NEXT J 
4255 NFYT 1 
4360 PRINT “FA 100 100 PU” 
1765 PRINTI "PA A -2 ER 1.5 2° 
4370. [tL GTCX min) 
4375 FAR Tel 10:51 
4780 Ya-Xmine10“(1-1) 
4386 X=LGT(Xa/Xmin)eSfFx 
4740 PRINT "PA”ıX .".0:” 
4338 IF 1i>=0 THEN PRINT “CP -2, 
4400 IF 14:0 THEN PRINT “CP -2.- 
440S 11=[1፥፦!| 
441@ NEXT I 
4415 PRINT "PU PA 0 0“ 
4420 Ti=LG6T(Ymin) 
4425  Y10-10 
4330 FMR T=1 TO (y+! 
4435 Ya=Ymin+10 (1-1) 
4440 ሃ።[ቦፐ(ሃጋ/ሃጠ1በ)*ፍያነሃ 
4445 PRINT “PA 3, 1.” 
4450 PRINT “CP -4,-.25:1.BIA;FR -2,2ILB"1T1447" 
4455 Tiselisl 
4450 NEXT 1 
4465 PFEP 
4370 INPUT “WANT USE 02۲۵۸۱۷۲۲ ٭ ]ع۱۵۸‎ (1=Y @=N)?* [dl 
4475 IF I414”>1 THEN 


SIP )ٴ۶‎ ۶۷۷۷۷۷۶۷٦ 
ILBIZSLER 8,26 1 ٭۳‎ 


4180  RFEP 
4485 INPUT “ESTER X-LAREL” XLahels 
4490 BFFP 


4495 TMHPUT “ENTER Y-LARFL" Ylabels 

4500 END IF 

4595 IF Ono 2 TIIFN 

451A PRINT "SR 1,.2:PU PA 40 -l4:” 

ASS ERIN 7 TER “5 3 İİ PRO 2.0 PUFPR 75,741: Buo1PR. 5.19” 
4820 PRINT “LR-TITR .5 -riLReatiPR .5.ta” 
46525 İF Tunza ዝክ 

4830 PRINT "LB)Y (K)” 

4635 1.2Е 

4640 FPIMT “IB) (F Y” 

46545 FND IF 

46550 ENN IF 

4555 (IF fipo=2 THEM 

4550 IF I[on=@ THFN 


ASRA PRINT "SR 1.5,2:PU PA 40. 154118 (Um SR | 1.Séi PR 0@.S 1]118215ፐ 1.ካ ۶۴ 
ہپ“ ۰ھ‎ 77۶ 
4510 ELSE 


E PRINT CSA | GS 2E EA 24-14 4L Bg CBtu/ZhriFR .5,.5:EB. :PR. e, En" 
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S ,21PR 7۰7 


46580 PRINT "LBftiPR .S,I1S5R | |,51[8ረ158 | 

4585 END IF 

4580 END IF = 

4595 IF 0no-0 THEN 

4500 IF IunsQ THEN 

4605 PRINT “SR |.S 2íPU PA -12 d@DI "ሾ,!፤[ኮባ (ህ/ጠ(ዞ8ጸ - ፀ.5158 |,|.,5([182158 | 
RE GE 

4610 ELSE 

4615 PRİNT “5R 1.5,2:PU PA -12.32:01 0,l1tBo (Btu/hriPR -.5 .SiLB.1PR .5,.517 
4620 PRINT ^LBft1S9R 1.1.5i1PR -1,.5sLB21PR. 1,.515R 1.5,21L B8)" 

4625 END IF 

4630 END IF 

4635 IF Qpo 0 THEM 

4640 IF 1ህበ=ፀ THEN 

4645 PRINT “SR 1.5,2:PH PA -12,3B:0I ወ.,|([(8ክ (W/MiPR -1,.SiSR 1,1.51LB2159R 1, 
5.2 .5,.51^ 

36504 PRINT '"[8.1ዞዞ .5.0([ርዞ'” 

4555 ELSE 

4660 PRINT “SR 1.5,2:PU FA -12:28:01 @,liLBh (Btu/hriFR -.5,.5;:LB. 1PR .5,.51° 
4665 PRINT "LBft14PR -1,.519R 1,1.81LB215R 1.5.21PR .S .S|6rLB. PR .5 ۳ 
4570 ENO IF 

4675 END IF 

4620 IF 14129 THEN 

4685 PRINT “SR 1.S.2 (PU PA S5@,-I6 СР" 1-1 ЕН ХТаһбе198)/2170018" 1Х1аһе|51”" 
4690 PRINI "ዞል -14 5ይ ርዞ 0.":-LENCY Iabe1$)0/2*5/541"0I ٥,۱۱۱۷۴" ۱ ۱۵ ٠ 
4695 PRINT "CP 60.0 DI“ 

4700 EMD IF 

4705 160 

4710 Reneat:1 

47157 ና]1=።] ፣+65 

4720 Yul--l.E96 

472 Icn-Q 

47:0 ВЕЕР 

4725 INPUT "WANT TO PLOT DATA FROM A FILE (1=Y O-N)?” O! 

4740 IF Ok=] THEN 

4745 BEER 

4760  IMPUT "ENTER IHE NAME OF IHE DATA FILE” O_files 

4755 ASSIGN @File TO O_file? 

4760 BECP 

326 ۶۳ 

477 INPUT “ENTER TIE BFGIHHING RUN NUMBER” Md 

37: REEF 

4780 INPUT “ENTER TINE MUMBEFR OF X-Y PAIRS STORED" ,Npairs 

4785 REEF 

4798 INPUT “CONNECT DATA WITH LINE (1-Y,.0-N)2" , Ic 

4:85 ٣ 

4800 PRIMTER 15 1 

4805 PRINT USING "4X,^"*Select a symbol:"^^ 

4BIO PRINT USING "6X,"*] Star ә Plus ion. 

4815 FRINT USING "6X "^3 Circle 4 Square” ”*” 

4820 PRINT USING "6X,"^S Romhus""" 

4825 PRINT USING "6X,"*65 Riqht-siıde-up trienale"”” 

4830 PRINT USING "8X,"^7 Up-side-doun triangle*”” 

4835 INPUT Svm 

4840 PRINTER 15 705 

4845 FRINT “PU DI” 

1850 IF Symzl THEN PRINT "SH*"* 

4855 IF Svms2 TIEN PRINT "SMe? 

4650 IF Svms3 THEN PRINT “Mo” 

486S IF Md: t THEN 

4870 FOR I=! 10 (Md-1) 

4875 ENTER $Filezira Xa 

4B80  HÉXT I 

4005 FEND IF 

4840 FOR Tel Tü Moaırs 

4895 ENTER BFileiYa Xa 
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4900 
4905 
4912 
4915 
4920 
“15 
4520 
4925 
4940 
4945 
4950 
4955 
4360 
4965 
4370 
4975 
4980 
4985 


4930 . 


4995 
5000 
5225 
5010 
5215 
5020 
5025 
59220 
5035 
5040 
5045 
5050 
5255 
care 
5055 
5070! 
5875 
5250 
5095 
5090 
5095 
5100 
5105 
5110 
SS 
5120 
ələ 
5130 
2136 
5140 
5145 
5150 
5155 
5166 
5185 
5170 
5175 
5190 
کے‎ 
5190 
25 
5200 
5:056 
5210 
Eos 
Soc 
S225 


THEN Ya=Ya/Xa 
THEN 


IF Opo=l 
IF Ünoe2 
Q-Ya 
Ya=Ya/Xa 
Xa=0 

END IF 
XacX11 TEN XlTeXa 
Xa»XuI THEN Xul=Xa 
Tun=? THEN 

ОосоД*2 THEN Xa=Xael.8 
Doo:0 THEN YasYas.176l 
069*8 THEN YasYas.317 
IF 00022 THEN Xa=Xae.317 
END IF 
Xx=LGT(Xa/Xminm)eSf. 
۷۰۱6۲۸۰۰ ۹۱۷ 
IF Svm-3 THEN PRINT 
IF Gen: A THEN PRINT 
IF IcI-0 THEN 

PRINT PA x v "° 
ELSE 
PRINT 
ENO ik 
IF $vm 3 
IF Svme4 
IF Svme5 
IF Sym=5 
If Syme? 
NEXT I 
PRINI 
BEER 
ASSIGN PFile TO a 
XII«*XII/!. 
XulsXuI*l. 
GOTO 8040 
END IF 
FRINT 
BEEP 
INFUT 
IF Go ርክ=! 
BEEP 
PRINTER IS | 

FRINT '!SING "4X ““Sel=ct 
PRINT USIHG "6x,"*"Q 
PRINT USING “5× ,۰۱ 
PRINT USING "BX,"*2 
INPUT Ion 

PRINTER IS 705 


17ت 


PSR EINAT TAS 


BAC K7 BD 


THEN PRINT 
THEN FRINT 
THEN FRINT 
THEN PRINT 
THEN PRINT 


“550126” 
"UC2.4.99,0,-8,-4,0,0,8 4,017 
"UCZ 0,99,-3,-6,-3,5,3,.6,3,-84" 
"UCG.5.3,99,3,-8,-6,0.3,81" 
"UC0.-5.3,99,-3,8,6,0,-3,-81" 


"PU" 


ГӘ TA 


"PU SM" 


“WANT TO PLOT A POLYNOMIAL (1*Y,Q-N)?" ,Go on 
THEN 


line type: 
Solid Tine 
Nashed””” 
5 Conger tine’ = dash 


REEP 

INFUT "SELECT (CsLIN,1*sLOGUX , Y ) D^ , ITog 
Ior nts! 

CALL Poly 


FOR x--0 TO C« STEP C./200Q 
Xa=Xmint+ 18 "Y. 

IF Xa XT] OR XaoXul THEN S200 
Ien=len+] 


۲ 


{E 
DE 
IF 
IF 
IF 
DE 
IF 
IF 
He 


Гоп= | 
Inn=2 
lon=3 
lon=4 
Топ=5 
Idf=1 
Оос=0 
Сосә2 
Оро=2 


THEN 
THFN 
THEN 
THEN 
THEN 
THEN 
THEN 


AND [log=1 


Idf2Icn 
Inf=Ten 
lgf= len 
Idf=Icn 
Idf*Icn 
Pu=l 


ሸበር 
MOD 
np 
Mon 
MOD 


^"2— ኤ የጋ 
nm 


Ya=fNPolv(Xa) 
AND Ilog=0 THEN Ya=Xxa/FNFolvíYa) 


[HEH ra=FNPoly( Xa) 
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5230 
5235 
5230 
5245 
5250 
5255 
5250 
9265 
5270 
Sits 
5280 
5285 
5290 
5245 
5200 
5305 
5310 
515 
JE 

55 
Seo 

S225 
5340 
5345 
5350 
53555 
5360 


u uu uvslun 
ed Cd Cji C4 L4 UL 

m 
ی ت ی ت‎ e ل٦‎ 


ул 
L4 
из ف‎ DID -) -J 
ርበ 


5400 
5405 
C410 
5415 
5420 
5425 
5430 
5135 
5449 
5445 
5450 
5455 
5460 
5465 
5470 
5475 
5480 
5485 
65430 
5495 
65900 
5505 
“512 
5515 
6520 
5525 
5530 
5535 
5549 
5545 
5550 


5555 


IF Ooo=1 THEN YasFNPoe]«(Xa)/Xa 

IF Ya'Ymin THEN 5200 

IF lun=1 THEN 

IF Ooo-2 ТНЕМ Хазхаз1.В 

1Е Ооо20 THEN Ya=Ya*.1761 

IF Opo=@ THEN Ya=Ya+.317 

IF Opo=2 THEN Xa=Xae.317 

ENO IF 

YəLGT(Ya/Ymin)əSfv 

کر ہر رک ری مھ مر وت 

IF Y O THEN ۷8 

IF Y>100 THEN GOTO 5300 

IF Pu=0 THEN FRINT "FA" .X,¥, "PU" 

IF Pu=1 THEN PRINT "PAÀ^,X.Y ,"PU* 

NEXT Xe 

PRINT "PU* 

EMO IF 

BEEP 

INPUT “WANT TO QUI] (15Y,05ND" ,Iat 

IF lgt=1 THEN 5335 

GOTO 4715 

PRINT “FU PA 0 0 SPQ" 

€ UBENÜ 

OEF FNHsmooth(X,Bco 1sat) 

OIM A(S) BCS),C( 5) OCS) 

DATA .20528,.25322..2190418,.55222,.79909 1.00258 
DATA .74815,.72992., . 73189. 7125709110 07 
DATA .41092..17726..285142,.54806 ..81316 1.0845 
OAIA .71403,.729132,.776585,.596581,.655867,.61889 
REAO Ale) Rie) Cte) Die) 

IF Boo 6 THEN I-8op 

IF Boo=5 THEN I=4 

IF Bop=108 THFN T=S 

IF Isat-=1 THEN 

HzsEXPCACIO*BCIP*LOG« X20) 

ELSE 

Has*EXP(K(GCCOIDO*DCI?*LOGCX D) 

END IF 

RETURN Hs 

FHENO 

DEE FNPolvOX) 

COM /Colv” A(10,.10),C( 10) .8(4) .Nop,Iornt ,Opo,IIog 
Xl=x 

Polv=B8(0) 

FOR I=] 10 Noo 

IF Iloo=1 THEN Xlet066X) 

PolvszPoly*BCID*X1^I 

NEXT I 

IF Ilon=1 THEN Polv=ExP(Poly) 

RETURN Poly 

FMEMO 

cue Poly 

OIM RETOS STAD SV ፍ- 12) 20100) YvC lO) 
COM ፡ርፀ1ሃ/ ፅ(19,|!ፀ) ር(|ቦ) B(4) N.İlornt ,Üco,llog 
COM /X«vy/ Xpn(5),Yot 5) 

FOR 1:0 TÜ 4 

В(1)=0 

NEXT 1 

Seer 

INPUT “SELECT (Q=FILE .1: FEYBOARO ,2=PROGRAM)” Im 
Im= I m4 1 

BEEP 

[ከቦ እ “ENTER NUMBER OF x- Y PAIRS" .No 

İF Im=l THEN 

BEEP 

INPUT “ENTER OATA FILE Hot" D file 
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5560 
5565 
5572 
5575 
5580 
5585 
5596 
5595 
5806 
5665 
58 | 0 
5615 
5620 
5625 
S630 
5635 
5640 
5645 
5650 
SESS 
5686 
SBS 
5670 
5675 
5622 
58ዩ5 
5639 
S695 
5700 
5765 
5718 
5715 
5720 
5725 


m 


sar 
5730 


5725 
5719 
5745 
5759 
5755 
S760 
5785 
5779 
5775 
5780 
5785 
5790 
5795 
5889 
5825 
5816 
5815 
5820 
5825 
5830 
5825 
5240 
5045 
S850 
5255 
5860 
5065 
5876 
5875 
5800 
5895 


ВЕЕР 

INPUT "LIKE TO EXCLUDE DAIA PAIRS (l=Y A=)?" led 
IF Ied2! THEN 

BEER 

INPUT “ENTER NUMBER OF PAIRS TD BE EXCLUDED" ,Ipex 
END IF 

ASSIGN ®File TO D file$ 

ELSE 

REEP 

INPUT “WANT TO CREATE A DATA FILE (1=Y 0@0=N)?" Yes 
IF Yes=1 THEN 

BEER 

INPUT “GIVE A NAME FOR DATA FILE” ,O_file$ 
CREATE BOAT O_ files S 

ASSIGN @File TO 09 

END IF 

END IF 

BEEP 

INPUT “ENTER THE ORDER OF POLYNOMIAL” .N 
FOR I*@ TO Ne2 

Sv(I)=Q@ 

S, (I )=Q 

NEXT ! 

IF Ied2! AND Ime! THEN 

FOR I=1 TO Ipex 

ENTER 8۲11 Y 

NEXT 1 

END IF 

FOR Tel TÜ fo 

IF Ims]! !HFN 

IF Ono=2 THEN ENTER 64۴ [168٤,۹۰ 
IF ዐፀፅፅ፡5 THEN ENTEP @Filery,.Xx 
IF Ilog=! THEN 

Xtzx/Y 

X=LOG(X 5 

Y=LOG( Kt) 

END IF 

ENO IF 

IF Im=2 THEN 

BEEP 

INPUT "ENTER NEXT X-v FAIR“ X,Y 
IF ves=1 THEN OUTPUT @F,il=ix . Y 
END IF 

IF Im 3 THEN 

X.(1)eX 

YvCIOsY 

ELSE 

XeXo(1-1) 

Y=Yp( 1-1) 

END IF 

8(02)=ነ 

5» (ሮይ )መ6ሃ(ወ)፥ሃ 

SCI»X 

5.(1)=ӱк(1)4Х 

FOR lel TÜN 

R(J)=R(J-1 )= X 

Sv(J)=Sv(J)*#R( J) 

NEYT J 

FOR J=2 10 Nee 

5(ህ)=5(ህ-1)5»፪ 

7717770 J ISK) 

NEXT J 

HEXT 1 

IF Yes=! AND Im=2 THEN 


BEEP 
PPINT USING "12% 00,"°" X-Y¥ pairs were stored in file 
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"" HOA" IND О files 


5830 
5895 
52900 
5305 
5310 
Sols 
5920 
2925 
5920 
5355 
5940 
5345 
5950 
5955 
S960 
5955 
5970 
S975 
5389 
5995 
5330 
5995 
6020 
6005 
5010 
6015 
5020 
6026 
6020 
6035 


END IF 

٥۵ھ(‏ )ء 5 

FOR IsQ TO N 
ር(1)=5ሃ፤![) 

FOR J=Q TO N 

QUL 13 ET) 
NEXT J 

NEXT 1 

FOR 1-0 TÜ N-1 
CALL Divıdet1) 
CALL Subtract( 141) 
NFXT I 
8(M)=C(N)/AIN N) 
FOR 1-0 10 М-1 
Bz 1.51) ፲ነ 
FOR JsQ TO I 


CON NDIE TOSA NE A ۱۱۲۰7۱7۱6۲110) 


MEXT J 


፲‹ዘ-1>1 )።8(8-1-1)/ጸ(አ-|-1 .ዘ-1-1) 


MEXT MI 
IPRINTER IS 701 

IFPINI B( +) 
IPRINTER IS 105 

IF Tornt=0 THEN 

PRINT USING “127. “”EXPCHEMT 
FOR I=Q TO N 


68۲٢۲٢۹٢٣٢ 


۴۱۲۷۰۱۱۲۲۰۰۰ SK MD. TOE SSI ECT) 


NEXT I 
FRINT ~ `“ 
FRINT USING "^12X.""DAT^ POINT 


6 = 


6040 
6045 
5050 
6055 
6060 
6065 
6272 
6075 
6030 
6085 
6030 
6296 
6100 
6166 
5110 
6115 
6120 
6125 
6122 
6125 
8140 
6145 
6150 
65 
6160 
6165 
5170 
6176 
6190 
6185 
6190 
6195 
6200 
6295 
6210 


FOR Tel TO Mo 
Ycs8(0) 

FOR J=) TON 
ኘር=መሃርሮ፥።1)።ዛኣ(!) J] 
NEXT J 

D=<Yv(I) Ye 


PRINT USING *15X,30,4X .4(MD.5DE , 1X0" 5 L,X& CIO, Yy CD), Yc D 


NEXT 1 

EMD IF 

ASSIGN @File TO ° 
SUBF ND 

SUB Divide(M) 


COM /Colv/ AC10,10)?,CC102,Bt 4D ,N, Iprnt ,Opo ,.L1og 


FOR [*M TON 
ጳፀ።ጳ( 1] Hi 

FOR J«M TO N 

ACT 17)5ሰሳ(1 .J)^Ao 
NEYT J 

Ct 1[IsC( 1)/Ao 
NEXT I 

SUBE ND 

SUB Subtract(k) 


COM /tolv/^ ^UI0.,10),CCIO) BC4O NlprnU”Uno 5 


FOR I*K TO N 

FOR J*W-] TON 
ACLUL.J)5ACK-1,J)-ACT .J) 
NEXT J 

6)) 0 

[IF I 

SUBF HD 

CUB Flin 


607/661۰۰۸0 6)۶" ۰ , Coe 0164 


COM /X«vy/ Xx(S59).Y. C5) 
PRINTER IS 705 
BEE 


10H 


YCCALCULATED? DISCR 


5215 INPUT “SELECT (@2h/h@% same tube, 1l=h( HF )/h( sm)” ,Irt 
62: BEEP j 
6225 INFUT “WHICH Tsat (1=6.7.07-2.2)" Isat 
.. Xmin=0 
23 Xnas^10Q 
6210 Xstep=2 
6245 IF Irt=@ THEN 
6250 Ymin=0 
6255 Yma.el.4 
6260 Ysten=.2 
ከ265 ELSE 
5270 Ymin=0 
6275 Ymax=15 
6282 Ysten=5 
6285 ENO IF 
56230 BEEP 
6295 PRINT “ItesSPtsIP 2300 2700 8300 583601" 
6300 PRINT "SC 0,.100,0 !@@:TL 72 01" 
6305 Sf»-z100/( Xmax-Xmin) 
6310 Sfv=100/(Yma--Ymin) 
5315 PRINT "PU Q.Q PO" 
B320 FOR Xa-Xmin TO Xma« STEP Xsten 
6325 X=(Xa-Xmin)+Sfx 
6320 PRINT "PA"ix,".Qi XT1^ 
5325 NEXT Xa 
E340 PRINT “FA 100 .0:PU:" 
6345 PRINT “PU PA 2.0 PD” 
6350 FOR YasYmin TO Yma« STEP Ystep 
6355 Y=(Ya-Ymin)»Sfy 
E360 PRINT “FA Q.":Y, YT? 
53585 NEXT Ya 
6370 PRINT “PA 0,100 TL @ 2" 
6375 FOR Xa=Xmin 10 Xma- STEP Xstep 
6389 Xe(Xa-Xmın )"6f. 
6395 PRINT “FA":X,".100: YT” 
£230 NEXT Xa 
6335 PRINT “PA 100,100 PU PA 100,0 PO” 
6400 FOR YarYmin TO Yma- STEP Ysteo 
6405 Y=(Ya-Ymin)+Sfv 
6410 FRINT "PD PA 100," Y."YT" 
6415 NEXT Ya 
6429 PRINT "PA 100,100 FU” 
BRS PRINT RAD SR L.S ”ج‎ 
6470 FOR Xa=Xmin TO Xmas STEP Xstep 
6435 Xə(2a-Xmin )eSf- 
6410 FRİNT “PA“:X.".0:" 
55212 EE İP” UP € ٢6٣٦ 
6450 NEXT xa 
6455 PRINT “PU PA 02.0" 
6460 FOR YasYmin 10 Ymax STEP Ysten 
6465 IF ABS(Ya) 1.E-5 THEN ٢٢۹ 0۵ 
5470 Y=(Ya-Ymin)»Sfy 
6475 PRINT "ቦጳ ዐ,"1ሃ.”” 
BASO ERUR CP -4 -.251LB' t Yar ^" 
6465 NEXT Ya 
56490 XIabeI$2'Oi1I Percent" 
6495 IF Irt=0 THEN 
5500 Ylabel$-"h/h0?7" 


5505 ELSE 

5510 Ylabel$="h/hsmooth” 

6515 END IF 

5520 PRINT "SR 1.5,2;PU FÀ 50 -10 CP":-LEHt4IabeI$)/21"0:L B^sXlabel$4"" 
Baz FRIHT "PA -11,58 CP 0,” :“(EHMYlebel$)3/2:5/6:”DI 0 .1:18”iYlabels$i:”” 
0520 PRINT “CP 2 0" 

6535 1ог=0 

6540 ВЕСР 
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8545 
5559 


5555 


6550 
6565 
5570 
6575 
6580 
6585 
6530 
6535 
6520 
6505 
6510 
6615 
5520 
6625 
6550 
6535 
6540 
6645 
6650 
6655 
6550 
6655 
6670 
5575 
6680 
6585 
6699 
6695 
6700 
6705 
5710 
5715 
6720 
67:5 
6720 
5735 
6740 
6745 
6750 
6755 
6762 
6765 
6770 

yıb 
5780 
6785 
6790 
5755 
5800 
6805 
5810 
6815 
6820 
6825 
6820 
6855 
6840 
6845 
6250 
6855 
6860 
6665 
6870 


INPUT “WANT TO PLOT DATA FROM A FILE (1=Y 0Q=N>?* Oko 
106-0 

IF Obo>=1 THEN 

BEEP 

INPUT "ENTER THE NAME OF THE DATA FİLE” .0.files 
BEEP 

INPUT "SELECT (Q@=LIMEAR , 1=LOG(X Y)° ,Ilog 
ASSIGN GFile TO D filef 

BEEP 

INPUT “ENTER THE BEGINNING RUN NUMBER” Md 
BEEP 

INPUT “ENTER THE NUMBER OF X-Y PAIRS STORED" ,Noairs 
BEEP 

INFUT "ENTER DESIREO HEAT FLUX",0 

BEEP 

PRINTER IS ! 

PRINT USING *4X,"*Select à svmbol:""" 

PRINT USING “4X.””1 Star 2 Plus sign... 
PRINT USING "4X ,""3 Circle 4 Gouarer?" 
PRINT USING "4X,^"5 Rombus""" 

PRINT USING "4X, "B Right-side-up triangle 5 
PRINT USING "44,777 Un-side-down triangle””" 
INPUT Sym 

PRINTER 1S 705 

PRINT "PU OI" 

IF Svms1 THEN PRINT "Che? 

IF Svm-»2 THEN PRINT *SM*" 

IF Svm=3 THEN PRINT “SMo” 

Nn*4 

IF IIogs! THEM Nn»! 

IF Md?! THEN 

FOR I=ft TO (Md-1) 

ENTER €Fiie1rXa Ya 

NEYT I 

ENO IF 

0150 

[ይ 1169"| [|ጀከ 4=!ዐ06(0ባ) 

FOR I=! TO Noairs 

ENTER @Fileska 8(፡) 

Ya=B(0) 

FOR Kk*1 TO Nn 

YasYa*B( b )+ O K 

NEXT K 

IF Tlogəel THEN YasENP(Ya) 

IF 1log-=9 THEN Ya=01l/Ya 

IF Irt=@ THEN 

IF Xa=0 THEN 

Yo=Ya 

Ya=] 

ELSE 

YaaYa/Yo 

END IF 

ELSE 

Hsm=FNHsmooth(0,Xa,Isat) 

Ya=Ya/Hsm 

ENO IF 

X.(I-1)=Xa 

ሃሃ(1-!| )=፥፣ሃ5 

Xs(Xa-*min)*Sf« 

Y*s(Ya-Ymin)*Sfvy 

IF Svm: 3 THEN PRINT "SM" 

IF Sym 4 THEN FRINT “SR 1,425 

PRİNT “PA X,Y °° 

IF Sym 3 THEN FRINT "SR 1.2,1.6" 

IF Svm=4 THEN FRINT “[/ር2,4.,99 G -8 -4.0 0 8 4.0 + ` 
IF Sym=S JHEN PRINT “UC3 Q 99.-3.-6.-3.5.3 .5.3.-61 
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5875 
5880 
6885 
6339 
6895 
6900 
6305 
5310 
591 
6920 
5325 
6320 
6925 
6940 
5345 
6350 
5 355 
6960 
6 365 
5370 
6978 
6380 
5985 
5390 
و‎ 
7002 
7005 
7010 
7015 
70-0 
"OS 
7626 
7935 
7040 
7035 
7050 
7955 
7060 
7055 
7076 
7075 
7680 
7885 
7030 
7095 
7100 
7105 
7110 
11) 
7120 
7125 
7130 
1155 
7140 
7145 
7150 
155 
7150 
7155 
7170 
roi 
7180 
7185 
7130 
7138 
7200 


۷۶۰58 < HEN ERINT UCI ۰۹9۰36 8" 
IF Syae7 THEN PRINT “UCU --5.3,33.-3.8.5.0.-3.-8:” 
NEXT 1 

BEEP 

ASSIGN @File TO s 

EMO IF 

PRINT "PU SM* 

BEEP 

INPUT “WANT TO PLOT A POLYNOMIAL ( 1e-Y .0-N)3?” lp 
IF Oko=f THEN 

BEER 

INPUT “SELECT (Q*LINEAR ,1-L0G(X ,Y))" ,Ilog 

lornt=1 

CALL Poiv 

FOR Xa-Xmin TO Kma, STEP Хз4ео/25 

Icn=Icn+1 

Ya=FNPol' (xa) 
Y=(Ya-Ymin)*Sfy 

Xc( Xa-Xmin)*Sf« 

IF Y Q THEM Y=9 

IF ¥ 100 THEN GOFO 7025 
Pu=Y 

IF lon=1 THEN Idf=Icnm MOO 2 
IF lon=2 IHEN Idf*Icn MOO 4 
IF Ion=3 THEN Idf=fen MOO 8 
IF Ton=4 THEN Idf=Icn MOD | 
IF TonəS THEN Idf=Icn MOO 3 
IF Idf*! THEN Pus! 

IF Pu40 THEM PRİNT "PA",X Y,"PO* 

If Puri THEM PRINT "PA" X,Y,"PU* 

NEXT Xa 

PRINT "PU* 

Ipnelontl 

COTO 6540 

ENO IF 

BEEP 

INPUT “WANT TO QUIT (1=Y дәнә?" Touit 
IF Iquit?! THEN 7070 

6010 6540 

PRINT “PU SPO" 

SUDENO 

SUB Stats 

PRINTER IS 701 

dağ 

K=0 

BEEP 

INPUT “PLOT FILE TO ANALYZE?" Files 
ASSIGN @File fO files 

BEEP 

INFUF "LAST RUN No?(Q0- QUIT)" .Nn 

IF Hn=0 THEN 7305 

Nn=Nn-J 

ӱк=0 

5 3 

$270 

S.5 =0 

Svs=0 

5:5=0 

FOR Tel TO Nn 

J=J+] 

آ, ء8۴۱۱1 ENTER‏ 

H-Q/T 

S. = 9. +0 
S-s=S.s +Q ` 
5ሃ*ሜዓኑ፥ፐ 
3vs=Sys+T 2 


rJ) o 


ዮጋ 
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Sz=Sz+H 

Szs=Szs+H*2 

NEXT I 

Оаҹе=5х/Мһ 

Tave=Sy/Nn 

Наҹе=52/Мп 

Sdevg=SOR(ABS 1 (Mn+*S-s-Sx 2 )/(ክአክበሰ*(ከበ-1)))) 
Sdevt*SQRCABSCUNn*Svs- Sy 2)/(Nns (Nn 12))) 
Sdevh=SOR(ABS1(Nn*S75-S2 2)/(Nn+(Mn-1))))> 
Sh=100+*Sdevh;Have 

Sa=100*Sdev0/Dave 

St=100«Sdevt/Tave 

IF K= THEN 7295 


PRINT 

PRINT USING "*1IX,""DATA FILE:^",1!4A" iF11e$ 

FRINT 

FRINT USING "11X,""RI'N Htube Sde vH Qdo Sdevğ Thetab SdevT”” 


ኮ= | 
FRINT USING "TIX DDy2 SX D.SDE Ix 30.20) 2X DOM. SD IK SO ED I Have SM Oe 


Go Taxe, BI 


7 300 
7265 
7310 
7315 
7320 
7329 
tose 
7355 
7240 
7345 
7350 
7255 
7360 
7255 
7276 
mir 
7360 
7385 
7330 
7,395 


6010 7115 

nSSIGCN 6ዩ11=ዱ1 [0 ፡ 

PRINTER IS 1 

SUBEND 

SUB Coef 

COM /Colv/ A(10,10).C( 102,Bt 4? ,N, Iprnt ,Qno, flog 
BEER 

INPUT “GIVE A NAME FOR CROSS-PLDT FILE” ,Cof$ 
CREATE BOAT Cof$,2 

ASSIGN @File TQ Cof% 

BEEP 

INPUT "SELECT «(Q*L INEAR,1-tL0GCX Y?)" [109 
PEEP 

INPUT "ENTER OIL PERCENT ¢-1=STOP)" 0 

IF Bep: 0 THEN 7390 

CALL Pol, 

QUTPUT 8Fi1lei:1Bop, Bt *) 

GOTO 7360 

ASSIGN 9፻11= [0 ə 

SUBEND 
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The 


Octagon 


APPENDIX C 
COHPUTER-CONTROLLED VALVE PROGRAM 


computer-controlled valve vas controlled by an 


Systems SYS-2A microcontroller. The S¥S-2A is a 


complete computer system on a single card and requires only 


a 5 volt supply for operation. The SYS-2A has provisions for 


4, 


2.9 


volt A/D data inputs and 8 high current digital 


outputs. The card was configured as shwon in Table 4 below. 


TABLE ü 
Configuration of SYS-2A I/O Channels 


Data rs Channel Purpose 
2#310 = 8 Valve position 
@#A10 = 9 Liquid R-114 Temperature 
Data E E Channel Purpose 
2#ል0 Shut Valve Relay (0=stop, l=shut 
ወ#ልዐ1 Open Valve Relay (0=stop, l=open 


The  proportional-integral-derivative control program 


below was written in NSC "tiny BASIC." The SYS-2A microcont- 


roller operator's manual includes an appendix of tiny EAS 


commands. Remarks are provided below that were not int 


IC 
he 


actual program (to minimize execution time) for clarity. The 


program was stored in EEPROM added to the card after program 
debugging was complete. 


10 
15 
20 
30 
40 
50 
60 


Let A=15:Let B=5:Let C=10:Let F=0:Let I=0 Initialize 
Let @#A00=O0:sLet Ә84А01=0 Stop Valve 
Print"Set Constants (1=Y,0=N) ?":Input Q 

If (0-0) then GOTO 60 

Berne 05 ,C,632":1nput HQU/E,L 

Let A=H:Let B=J:Let C=K: Let G=L 

Print"input Desired Temp(-1=End) ":Input R Enter Temp 
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70 Let @#A10=9:Delay 2:Let M=16*0#A11+0#A12/16 Read Channel 9 
80 Let @#A10=8:Delay 2:Let P=16*O#A11+0#A12/16 Read Channel 8 
90 Print"Tenp=",4,"Pos=",P 


100 If (R<0) then GOTO 800 Exit Program 
200 Let @#A10=8:Delay 2:Let P=16*0#A11+0#A12/16 Read Channel 8 
210 Let S=0 Ave 10 Readings 


220 For N=1 to 10 Step 1 
230 Let @#A10=9:Delay 2:Let M=16*@#A11+#+0#A12/16 Read Temp 


240 Let E=M-R: Let S=S+E Error; average 
250 Next N 

260 Let E=S:Let D=S-F:let F=E Error; Derivative 
262 If(E<0)then I=I-1 Integral 

264 If(E>0)then I=I+1 Integral 

266 If (I<-31000) then I=-31000 Binary 0۲ 
268 1f(1531000) then 1-231000 Binary Lime 
270 PrintuiM=W ا ا‎ ٤> ELO p MISS p 

300 If (E<-1000) then GOTO 600 Control Band 
310 If(E>2000) then GOTO 700 Control Band 
400 Let V=(E/A)+ (B*D)+ (1/C) PID Value 

410 FrintwWSy="wu, v ,"w"gT=" E k "pDISwW B*#p ə 6 Debug Tool 
820 If(V»0)then GOTO 500 Above Desired 
430 Let V=V* (-1) *G Shut Faster 
440 If(V<1) then V=1 Check 


450 If(V>1040) then V=0 Delay Function limited 1-1040 entry 
460 If(P<150)then V=15 Slow Valve near end of travel 

465 If(P<50)then =10 

470 If(P<20)then GOTO 200 End of Travel Limit 

480 Let @#A00=1:Delay(¥):Let @#A00=0 Shut Variable Time 
490 GOTO 200 

500 If(P<2500)then GOTO 200 End of Travel Limit 


510 If(V<1)then V=1 Check 
520 If(V>1040) then V=0 Delay Function limited 1-1040 entry 
530 Let @#A01=1:Delay(V):Let @#A01=0 Open Variable Time 


540 GOTO 200 
600 If(P€20)then GOTC 200 End of Travel Limit 
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610 
620 
630 
640 
645 
650 
660 
670 
700 
710 
720 
730 
790 
800 


Eriınt 600 Loop" Out of Bound Warning 


Let V-0 Max valve shutting speed 
If(P<250)then V=50 Slow near end of travel limit 
If(P<150)then V=15 


If (P<50)then V=10 

Let @#A00=1:Delay(V¥):Let @#A00=0 Shut Valve Fully 
Let @#A10=8:Delay 2:Let P=16*@#A11+0#A12/16 Read Channel 8 
GOTO 600 

If(P52500) then GOTO 200 End of Travel Limit 

Print"700 Loop" Out cf Bound Warning 

Let @#A01=1:Delay O:Let ወ#እዐገ=1 Open Valve Fully 
Let 3$A10-8:Delay 2:Let P=16*0+*A11+0+*A12/16 Read Channel 8 
GOTO 700 

END 


APPENDIX D 
EXAMPLES OF REPRESENTATIVE DATA RUNS 


The two printouts below 


are samples of the output of 


DRP2. The first printout (WH79) 


high heat flux with 0 percent 


is £or the smooth tube at 


OL The second grintout 


(HF101) is for the High Flux tube under similiar conditions. 
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Month, date and time :02:25:10:02: 
NOTE: Program name : ORP2 
Disk number = 727 
Old file name: WH79 
This data set taken on : Q2:10:12:50:15 
Tube Tvoe: 4 
Data Set Number = 1 Bult Oil X = 0.0 02:10:13:16:18 
TC No: 1 2 3 4 5 5 7 8 
Temo : 16.42 14.69 14.61 15.00 14.67 14.91 16.79 14.74 
Tua Tliod ۲1۱:892 Tvaor Psat Tsumo 
15) 92.25 -2.120 -.63 -6.31 -17.?7 
Thetab Htube Odo 
16.935 5.704۲09 3 64 
Qata Set Number * 2 Bulk O11 X = 9.9 ص2۰۶5‎ 
TC No: | 2 3 4 S 6 7 8 
Temo : 16.48 14.60 14.54 14.33 14.52 14.74 16.69 14.87 
Tura Tliad Tliod2 Tvaor Psat Tsumo 
15415 “0 -2.10 -1.44 -6.32 -17.1 
Thetab Htube Odo 
16.866 S.745E+03 3.589E+04 
Data Set Number = 3 Bulk 011 2%» 0.0 02:10:12:13:38 
TC No: l 8 3 4 2 5 7 3 
Temo 16.48 14.56 14.52 14.35 14.50 14,67 16.62 14.67 
Tua Tliod Tliad2 Tvaor Psat Tsumo 
15020”? -2.13 -1.51 -8.32 -17.0 
Thetao Htune Qdo 
16.840 5.760E+03 3.701E+04 


NOTE: 03 X-Y vairs 


DATA FILE:PuH73 


Odo 
9.683E+04 


SdevH 
.21 


RUN Htube 
35.7366٤03 
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Sdev0 


were stored 1n plot data file PUH79 


Thetab SdevT 
16.880 :29 | 


ማጫ 
٭‎ Ae سے‎ 


Month, date and time :03:<5:10:03:26 
NOTE: Program name : DRP2 

Disk number = 02 

Old file name: HF101 

This data set taken on : 92:16 


Tube Tvpe: 5 


Data Set Number = 1 Bulk Oil X = 


TC No: 1 2 3 4 
Temo 25 DU DS ə 4.57 
Twa Tliad Tliad2/ Tvapr Psat 
(a mm Zal 21287 ን BUR 
Thetab Htube Qdo 
3.370 2.801E404 9.4406 +4 
Data Set Number = 2 Bulk Oil X = 
TC No: | ረ 3 4 : 
Tema 5.54 7:5: 3.76 4157 
Tua Tiiad Tliqd2 Tvaor Psat 
15575: 4 -1.90 22: “5:0 
Thetab Htube Qda 
2377092926+ 049۳4436094 
Data Set Number = 3 Bulk Oil X = 
TC No: | 2 3 4 
Tema 555 28: 3.75 4.57 
Twa Tliad Tlıad2 Tvaor Psat 
4256.02.14 =) "ር 9... 
Thetab Htube Oda 
5195121115) ሃለወ 4 ደ 45БЕ 104 


513:2552 
0.0 +0 + + +۶ 
5 6 7 8 
3.93 4.52 7. 31 4 17 
Tsuma 
-15.7 
0.0 02:15:19:31:15 
5 5 7 8 
2:05 4.53 T3] 4.18 
Tsump 
-15.6 
0.0 de: 16:19:31:27 
5 6 7 8 
3.93 4.53 ə 42:17 
Tsumo 
-15.5 


NOTE: Q3 X-Y pairs uere stored in plot data file PHF101 


DATA FICE: PRETO) 


Qda 
9.443E+04 


SdevH 
ər 


RUN Htube 
3 94 
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SdevQ 


Thetab SdevT 
25 30755 .15 


. | ә 


APPENDIX E 
UNCERTAINTY ANALYSIS 


The uncertainty of the heat-transfer coefficient at 3? 
kW/m? and 5 kW/m2 of runs WH79 and HF101 are analyzed below. 
The analysis is based on the Kline-&McClintock [Berle 
method of uncertainty analysis. 


The heat-transfer coefficient is: 


hmc dco ни (E.1) 
Tuo ^ Tsat 2 
and | 
DA ə Twi - m - feat (555 
where 
h » heat-transfer coefficient 
gar heat flux corrected for end losses 
Two > average outer wall temperature 
Tsat = Saturation temperature 
ru - average inner wall temperature 
2 heat input corrected for end losses 
D2 = outer wall diameter of tube 
Dj = inner wall diameter of tube 
k = thermal conductivity of wall 
L = length of heated surface 
let 


- ዐር ይቪ (D5/D1) 
: ни (Е.3) 
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According to Kline and HcClintock, the uncertainty in the 


kaoat transfer Coefficient is: 











2 _ 2 
h qc Tuo ^ Tsat 
2 27:12 
a OF E ÖTsat 
Tuo 7 ፲5ኋ፻ Two - Tsat 


by neglecting the erter from the Logarithmic term, because 
1t is small compared to the other terms, the uncertainty of 


the Fourier term (F) can be estimated as: 








" 2 2 2 1/2 
“EF 57:17 HE (8.5) 
Qc k E 
and 
Oc = qe TT 1) L, (E. 6) 
which has an uncertainty of 
2 2 2 112 
80 6 8D 
-. | 4(-2 s (E .7) 
e ባር 03 L 


Table 5 lists the various terms of equations (E.1) through 
(E-7) assuming all 8 wall thermocouples are used to get the 
average wall temperature 135 The data-reduction prcgran 
(DRP2) used this method to calculate the heat-transfer 


LT 


TABLES 


Uncertainty Analysis Terms Using 8 Thermocouples 





File ИН79 4H79 HF 101 HF101 
Heat Flux 37 kW/m? 5 kw/m? 37 kW"/m? 6 kV/n 2 
م۵‎ 
nr 0.007 0.009 0.006 02/00) 
Sk 
ES 0715 OS 0:33 02-39 
SL 
— 0.0005 0.0005 O50005 0. 0005 
F term (°C) 0.201 0.024 150106 0. 170 
SF (°C) 0:030 0.004 0.434 0.056 
qu OC) 11.46 53 0.79 -1.43 
vi 
T (ОС) =2.21 -2.23 -2.21 -2. 18 
Sat 
Tu; 
—— 0.032 0.073 0.290 0.12: 
Too - Tsat 
5 0.002 0.0005 0. 259 0.097 
m B Tsat 
“€ 0.003 0.008 0.015 0.019 
Two 7 Teat 
Öqc 
: 0.003 0.007 0.008 0.005 
e 
Oh 
S 0 2 O07 4 0.389 0. 160 
h (W/m2 K) 2660 540 22300 8400 


coefficient. Table 6 lists the same terms assuming only the 
center 4 wall thermocouples (2, 3, 25, and 6) are used to 
calculate the heat-transfer coefficient. Comparing Tables 5 
and 6 shows the effect of the axial and circumferential wall 
temperature distributions on the uncertainty of the 
heat-transfer coefficient. 

The constraining error of the smooth tube is the uncer- 
tainty in the wall temperature. Removing the effect of the 
axial wall temperatue distribution reduces the uncertainty 
by 1 percent. The large wall superheat of the smooth tube 
contributes to the small magnitude of the uncertainty terms. 
The constraining error of the High Flux tube is also the 
uncertainty in the wall temperature, but the uncertainty of 
the thermal conductivity (part of the F term) is about the 
same magnitude an results in the larger overall uncertainty 
of the High Flux data. The small wall superheats of the 
High Flux tube also amplify the magnitudes of uncertainty 
terms. Removing the effect of the axial wall temperature 
distribution makes the uncertainty in the thermal conduc- 
tivity of the copper-nickel High Flux tube the constraining 
uncertainty. The axial temperature distribution is respon- 
Sible for about 25 percent of tke wall uncertainty term, but 
again the combined effect of the uncertainty of the thermal 
conductivity, and low wall superheats, does not make it 
wholely responsible for the large uncertainty of the High 
Flux data. More accurate data could be obtained on the High 
Flux surface by using a solid copper tube without a large 
uncertainty in the wall resistance or an axial wall 
temperature variation. 
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TABLE 6 


Uncertainty Analysis Terms Using Center 4 Thermocouples 








File WH79 WH79 HF101 HF101 
Heat Flux 37 kW/m? 5 kWi/n? 37 kW/m? 5 kW/m? 
60, 
Bor 0.007 0.009 0.006 0.007 
дур 
E 0.15 0.15 0.33 0.33 
SL 
-. 0.0005 0.0005 0.0005 0.0005 
F term (°C) 0.201 0.024 1 6 o 
5Р (°C) 0.030 0.004 0.434 0.056 
T (9C) 11.34 5.80 0.44 -1.49 
WI 
T (9C) 22.24 0005 2 E 
Sat 
6T. 
a. 0.015 0.063 02219 0.097 
Tuo ^ Tsat 
SF 
0.002 0.0004 0.325 0.108 
Tuo = Tsat 
Tsat 0.003 0.008 0.019 0.021 
Tuo с Tsat 
5 
de 0.003 0.007 0.008 0.003 
qc 
дә 
ES 0.016 0.064 0.389 0.146 
h (W/m2 K) 2690 530 27800 9280 
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